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FOREWORD 


Efficient management of the Large Area Crop Inventory 
Experiment (LACIE) dictates that effective controls of proj- 
ect activities be established. To provide a basis for effec- 
tive control/ documentation will be prepared/ baselines will 
be established/ and changes to the baseline will be subse- 
quently controlled by the proper management levels. 

The specific control documents which will be used are 
defined in the LACIE Project Plan/ Document Number LACIE- 
C00605. All elements of the LACIE project must adhere to 
these baselined control documents/ and/ where it is consid- 
ered that the requirements should be changed/ the proper 
change request accompanied by a full justification must be 
submitted to the proper management level in accordance with 
established procedures. These documents will be maintained 
current by change notices and revisions as required. Each 
change notice and/or revision will reference the applicable 
Change Control Board Directive {CCBD) which approved the 
change . 

This document/ LACIE-C00200 / Volume II/ Revision E, de- 
fines the LACIE Classification and Mensuration Subsystem 
(CAMS) requirements and has been prepared in accordance with 
the "Instructions for Preparation of LACIE Requirements 
Documents"/ LACIE-00100/ Revision C, dated November 20/ 1974. 
"Full-Up System" as used in this- document is defined as the 
system required to accomplish LACIE Phase II. In general/ 
the approach used in each section is to first specify the re- 
quirements of the "Full-Up System" and then to specify the 
requirements of any interim systems by reference to specific 
paragraphs in the Full-Up System requirements sections of the 
document. The LACIE Project Phases are defined in the LACIE 
Project Plan/ IA.CIE-C00605 . 

The organization responsible for the implementation of 
eaph requirement defined in this docviment is specified on an 
individual requirement basis. Where the implementation re- 
sponsibility applies to the complete section/ the implementa- 
tion responsibility is specified after the section title. A 
"section" for the purposes of designating implementation 
responsibility is defined as being any numbered paragraph 
and all subparagraphs. Where different implementation re- 
sponsibilities apply to different portions of a section/ the 
implementation responsibility is specified on an individual 
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paragraph' or sentence basis, as applicable. All implementing 
organizations designated shall accomplish their implementa- 
tion activities in accordance with the requirements specified 
herein . 




R, B. MacDonald 

Manager, Large Area Crop Inventory Experiment 
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SECTION 1.0 

FUWCTIOUAL RESPONSIBILITIES 


1.1 


1.2 

1 . 2.1 

1 . 2.2 

1.2.3 

1.2. It 


GENERAL 

Design, implement, and operate a classification subsystem to de- 
termine values for the proportion of wheat in each Landsat acquired 
sample unit defined by the Crop Assessment System (CAS), at each of 
as many as foior periods in the growth of wheat. The subsystem 
shall be designed to utilize Landsat data as the primary source. 

SPECIFIC 

Subsystem Implementation 

Develop, specify, and implement the hardware and software systems 
required to meet the operational requirements of section 1.2.4. 

Personnel 

Identify and tram the personnel needed to operate the above 
system in accordance with the procedures of section 1.2.3. 

Procedures Development 

Develop procedures to accomplish the operational responsibilities 
of section 1.2.4 with the system of section 1.2.1 and the personnel 
-'of section 1.2.2. 

System Operations 

Estimate, for each acquisition of each sample unit, the proportion 
planted in wheat and, for those acquisitions occurring after wheat 
harvest, the proportion of wheat which was harvested. Estimate the 
error in each of the above estimates. Detennine the resources con- 
sumed in the above operations. Report all of the above to the CAS. 
Support CAMS Post Operation Performance Analysis by Testing and 
Evaluation. Perform additional analyses of sample segments in 
support of accuracy assessment and Test and Evaluation. 
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SECTION 2.0 
APPLICABLE DOCUMENTS 


The following documents are applicable to the extent specified 
within : 

A. Large Area Crop Inventory Experiment Project Plan - 
LACIE-COO 605 , August 1975 . 

B. Earth Resources Interactive Processing System Requirements 
Document PHO-TR514, June 197^ 

C. Flowcharts for LACIE Operational CAMS, December 2, 197^ 

D. LACIE CAMS Training Plan LACIE-OO60I, April 1975 
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, SECTION 3.0 


3.1 SUBSYSTEM DEVELOPMENT 

Figure 3-1 indicates the task flow required to place the CAMS into 
operation. 

3.2 SUBSYSTEM OPERATION 

Figure 3-2 shows the basic CAMS functions in terms of gross inputs 
and outputs. 

3.2.1 Plan 

Operation of the CAMS is based upon multispectral/multitemporal 
pattern recognition of Landsat data. In this process, the identity 
of the ground cover class for each Landsat picture element ("pixel") 
is estimated by comparing its radiance in four spectral bands for 
each of as many as four dates to the distributions of radiance 
(versus spectral band and date) for various candidate ground cover 
classes. In this classification process, the distributions are 
assumed to be multivariate normal, and the decision as to the most 
likely, class is made in a maximum likelihood fashion. 

Such a classifier requires statistics (mean vectors and covariance 
matrices) which describe the radiance distributions for the various 
candidate classes. In the CAl©, such statistics are estimated by 
computing the means and covariance matrices from actual Landsat 
data over earth surface areas covered by the actual ground cover 
classes of interest. 

These distributions ("signatures") are functions of many variables 
such as soil color and moisture content , growth stage of plants , 
atmospheric haze and water content, and solar elevation angle; 
therefore, in the CAMS, the statistics must be computed for regions 
in which these parameters are similar to those expected in the 
region to be classified. Additionally, the CAMS provides capa- 
bilities of correcting for changes in solar elevation angle and 
atmospheric effects. 

Two approaches will be used to find suitable areas in which to 
train the classifier prior to classifying a particular segment of 
Landsat data. In the first approach, analysis interpretation (AI) 
of imagery produced from Landsat segment data will be performed to 
find areas of wheat and non-wheat ground cover classes. Statistics 
will be computed over these areas ("training fields") and used to 
classiiy other pixels in the same Landsat segment or in other 
nearby segments. In the second approach, "analog regions" in 
which the signatures of wheat and non-wheat are highly similar 
to those in the recognition segment will be located. Then, sig- 
natures will be computed for areas of known ground cover classes 
and extended to the desired recognition area. These processes are 
indicated in figure 3-3 and 3-^. 
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iCAMS development flow. 
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Figure 3-t| , Illustration of analog site signature extension. 


Image analysis to locate training fields -will not be performed in 
all segments. It is more cost effective to locate training fields 
in an effective subset of the segments and extend the signatures 
to all segments . It is estimated that this subset of the segments 
("training segments") will contain 20 percent of all segments. 

Image analysis will also be used for another purpose. Previously, 
difficulty has been experienced in accurately classifying such non- 
cropland areas as forests, and towns. Frequently, such areas are 
highly inhomogeneous at the Landsat resolution, and poor classifi- 
cation has resulted. Therefore, in the CAMS, major noncropland 
areas will be identified by image analysis and will be manually 
excluded. Similarly, areas in which data quality is poor due to 
such factors as clouds and noise can be identified manually . 

Wheat-growing regions of the earth will be partitioned into 
smaller areas over which signature extension is expected to be 
successful. This partitioning will be accomplished based upon 
such ancillary data as crop calendars , meteorology , and soil 
color, as well as on the basis of trial classifications. Each 
signature extension region will contain one or more training seg- 
ments. If cloud cover or other operational problems prevent acqui- 
sition of the training segment for a given signatiare extension re- 
gion, signature extension will be attempted from a neighboring 
region. This approach is expected to provide reduced, but still 
acceptable , classification accuracies in the region without a 
training segnent. It is anticipated that signature extension re- 
gions typically will have approximately the area of a Landsat frame. 
This approach is also indicated in figure 3-3. 


The CAMS design is based upon the estimate that wheat can be sepa- 
rated adequately from other crops by analysis of up to four acquisi- 
tions of Landsat data during the growth of wheat . The Landsat ac- 
quisitions are to occur in a sequence which maximizes the discrimi- 
nation of wheat from non-wheat and the usability of the obtained 
inforBiat ion . 

Nominal Biowindow Definitions (Historical Data) 

A . Biowindow I 

1. Winter Wheat (South of N lat) - U.S., U.S.S.R. 

Opens one month past when 50 percent of the strata was 
planted or at stage 2.3 (when crop is 30 percent between 
emergence and jointing). In other wording, stage 2.3 is 
estimated as 30 days after 50 percent planting date. Close 
is one day before 50 percent of the strata has been reported 
in the jointing stage. 

2. Winter Wheat (North of 4U® N lat) - U.S., U.S.S.R. 

Open same as above. Closes when snow cover dominates the 
image or when climatological average one inch snow line 
occurs or about 50 days after stage 2.3. 
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B. Biowindow IB 

Winter Wheat (Worth of 44*^ If lat) - U.S., U.S.S.R. 

Opens 50 days Before 50 'percent jointing (based on nominal 
SRS reporting data for 50 percent jointing) . Close's one 
day before 50 percent of the strata has been reported in 
the jointing stage . 

C . Biowindow I 

Winter Wheat - Brazil, Argentina, Australia, India, China 
Opens at the 50’' percent planted date. Closes one day be- 
fore 50 percent of the strata has been reported in the 
jointing stage, 

D . Biowindow I 

Spring Wheat 

Opens 15 days before 50 percent of the strata has been 
reported in the planting stage. Closes one day before 50 
percent of the strata has been reported in the jointing 
stage . 

E . Biowindow II 

Winter Wheat and Spring Wheat 

Opens when 50 percent of the strata has been reported in 
the jointing stage. Closes one day before 50 percent of 
the strata area has been reported in the heading stage. 

F . Biowindow III 

Winter Wheat and Spring Wheat 

Opens when 50 percent of the strata has been reported in 
the heading stage. Closes one day before 50 percent of 
the strata has been reported in the soft dough or turning- 
yellow stage . 

0 . Biowindow IV 

Winter Wheat and Spring Wheat 

Opens when 50 percent of the strata has been reported in 
the soft dough or turning-yellow stage. Closes l8 days 
after 50 percent of the strata area has been reported 
harvested. 

For LACIE Phase III, the above also is acceptable for all training 
segments . 

For LACIE Phase III, window IB will be used only for those segments 
for which no satisfactory interpretations were made in window I. 

All segments will be acquired by GSFC at every available opportu- 
nity; however, to reduce the cost of preprocessing, transmitting, 
and storing the data, ordinary segments will be preprocessed and 
transmitted to JSC only once (if acquired) per stage. Designated 
training segments will be preprocessed and sent to JSC for each 
acquisition lying within a development stage. This is indicated 
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in figiire 3-5 in ‘whicli acquisitions of ordinary segments are 
indicated by X's and acquisitions of designated training segments 
are indicated by circled X's. The purpose of the extra acquisi- 
tions is explained below. 

Signatures obtained on one acquisition within a -development stage 
are not necessarily expected to be valid for other acquisitions 
within the same phase; however, training fields obtained on one 
acquisition will usually be valid on other acquisitions within the 
same phase. Therefore, by reacquiring the training segments and 
using training fields located on previous acquisitions, statistics 
appropriate to an acquisition can be computed without expending 
image analysis on that acquisition. It is planned to conduct major 
image analyses only upon training segments and only once per devel- 
opment stage. 

Fxgure 3-5 also indicates, in the column for segment 7> that where 
no suitable training data are acquired, it will be necessary to 
perform additional image analysis . Such analysis will include not 
only the most recent acquisition but also previous acquisitions. 
Changes in training segment selection will be made \mder procedural 
control . 

Selection of training segments for use in classifying individual 
recognition segments will be a dynamic process . Figure 3-6 indi- 
cates the inputs, considerations, and outputs in this process. 

3.2.2 Implementat ion 

Figure 3-7 gives a superficial view of the data processing required 
to implement the foregoing plan. Incoming Landsat data (which 
arrive at JSC one acquisition at a time, but are preregistered so 
that ground positions on different acquisitions bear the same scan- 
ner- coordinates) are used to generate film for image analysis and 
are input to the computer for composition and indexing (C&l). The 
C&I application stores the newly acquired data and prepares indexed 
data reports. (The C&I application is not discussed further in this 
volume, as it is a part of the Information Storage, Retrieval, and 
Refonnatting Subsystem (ISRRS).) Working with the film and the in- 
dexed data reports, analysts enter field coordinates and job con- 
trol input to the computer for classification processing. Completed 
classification jobs are screened for classification quality, and, 
where necessary, jobs are reworked. 

Figure 3-8 shows that the classification processing will be operated 
in two modes, a semiautomatic (batch) mode and an interactive mode. 
The major function of the data processing is expected to occur in 
the hatch mode, with the interactive mode being reserved for those 
jobs not practicable in the batch mode. 

Figure 3-9 shows a more detailed view of' the image generation 
process. Incoming Landsat data are composed with previously 
acquired Landsat data, allowing the generation of multitemporal 
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Figure 3-6.- [Training segment selection. 
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Figure 3-8.- [Functional overview of LACIE CAMS 
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ADP RESULTS 



Figure 3-9.- ,I>ACIE CAMS functional image flow 



images. Additionally, the composed data are returned to the image 
data base where they are available for use in classification proc- 
essing. Film is generated -using the Production Film Converter (PFC) 
in building 30, or as a backup, the Aerojet Data Analysis Station 
(das) in building 17. Film is transmitted to the JSC Photographic 
Technology Division (PTD) for development and then delivered to the 
Earth Observations Division (EOD) for analysis. Note that some ADP 
results are also transmitted in the form of imagery. 

Figure 3-10 shows the flow of field definitions obtained by Analyst 
Interpreters (AI) and Data Processing Analysts (DPA) . Field defi- 
nitions, in scanner coordinates, are input, via a work control 
center, to the fields data base in the ISRRS (formerly IMS). The 
ISRRS issues, via the PFC and the PTD (not shown), a transparency 
at the same scale as other LACIE imagery showing the boundaries 
and identities of all fields . 

Figure 3-11 is a more detailed -view of the batch mode classifier. 

The interactive classifier is similar, except that the analyst, 
working at a terminal, replaces the CAMS supervisor. In the classi- 
fier, training field coordinates, image data, and instructions are 
input to the statistics computation process. Following statistics 
computation, modifications to statistics to account for changes in 
siui angle or atmospheric turbidity can be performed. Computed 
statistics are available as an output at this point; alternatively, 
external statistics can be input at this point. This allows the 
development of offline statistics modification programs. The 
statistics can then optionally be input to the feature extraction 
processor in which those Landsat channels containing the most 
information can be extracted. Next, image data from the recog- 
nition segment and statistics are input to the classifier, which 
classifies each pixel in one of the input categories. The 
classifier offers the option of summing the density ftinctions of 
a number of' classes and/or subclasses into a category density 
function. 

Following classification, certain critical results are input to 
the ISRRS where they are accessible to the CAMS and other subsys- 
tems. More extensive results, including diagnostics and all data 
reguired to duplicate the run, are output to the analyst. Follow- 
ing evaluation of the resTolts, the analyst stores an evaluation of 
the results in the ISRRS. 
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Figure 3-10.- ILACIE CAMS field definition flow 
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Figure 3-1 1 LACIE classification 
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SECTION Ij.O 

CLASSIFICATION MD MENSUEATION SUBSYSTEM REQUIREMENTS 
FULL-UP SYSTEM 
Data Systeias Requirements 

(Req.’d "by CAMS; Cat. 1, except as noted; Impl. Resp. : CAMS/DSAD) 

This subsystem is enTisioned as several applications programs 
capable of being operated in several modes. The design is required 
to be modular so as to allow subsequent additions and changes in 
algorithms and procedures with minimum delay. 

The system shall be a production line system consisting of a rela- 
tively inflexible semiautomatic mode and a highly flexible inter- 
active mode, both of which are of operating on a common data base. 
This data base will contain extensive Landsat segment data. While 
it is not required that this entire imagery data base be simultane- 
ously accessible, it is required that some means be provided to en- 
sure that those segments to be processed together be simultaneously 
accessible when needed. 

The system must be designed so that the following conditions are 
met : 

A. An average throughput of (420 Phase I, 7l4 Phase II, and 840 
Phase III) sample segments/month can be processed, after 
allowance for downtime. 

B. Turnaround time on each processing job is less than 2k hours. 

C. Any segment can be processed on any day, regardless of the 
processing status of other segments . 

Status reports indicating the status of all Landsat segment data 
in the system are required. It shall be the responsibility of the 
EOD to verify performance of processing jobs . 

The semiautomatic mode should be designed to minimize manual 
handling, particularly transcription, of processing instructions. 

This Requirements Document references flow diagrams which are 
intended to clarify the verbal description of the processing. 

While these flow charts necessarily imply an implementation, the 
implementation shown should be considered only illustrative. Any 
other implementation which is functionally equivalent is clearly 
acceptable . 

Flexibility shall be provided to allow change of control parameters 
in the semiautomatic mode under configuration control. This appli- 
cation program must be able to perform classification on any number 
of categories up to 26 . Each category may be divided into one to 
, 60 classes. Each class may be subdivided into one to 60 supclasses. 
The total number of subclasses for all categories combined shall not 
exceed 60 . 



U. 1.1.1 Channel requirement s .- Purging and retention of these data .are TED, 
hut ciirrent planning should he based upon the assumptions that 
image data will he retained for active use three months after the 
end of a ‘crop year ‘and that purging and" reinitiation will occur 
upon a crop year basis. 

During this TED period, data from two crop years may exist in the 
active image data base for some segments.’ The capability to pm:*ge 
either year is required. This option must be available for 
Phase II. 

Channels and acquisition dates to be used in further processing 
must be selectable on a per segment , by user input . 

1;.1.1,2 Fundamental unit .- The fundamental unit of data to be processed at 
each distinct execution of this application program shall be a 
sample segment or a collection of sample segments to be classified, 
and a set of one to three training segments. Member (s) of the 
group of sample segment (s) to be classified may also be members of 
the group of training segments. 

Fields .- The user shall be able to identify the following types of 
fields in each segment : 

A. Class training fields (up to -60 subclasses) 

E. Class test fields (up to 60 subclasses) 

C. Designated Other (DO) fields 

D. Designated Unidentifiable (DU) fields 

E. Mean Level Adjustment (MLA) fields 

Data base : 

A. A fields data base, which may be vacant, shall exist for every 
segment , and contain up to l46 fields per each segment . 

B. A dot data base shall exist which will contain a set of fixed 
single pixel locations, with the category level designation and 
function designation for each segment . The data base is updated 
by an offline utility program. 

C. The capability should exist to store multiple sets of dots and 
fields for each segment so that any set of dots and the corre- 
sponding set of fields can be used to process a segment . 

■When a segment is represented in the active image data base by data 
from two crop years, separate fields data bases shall exist for the 
two separate years. 

These fields shall be used as follows : , 

A. Training fields will be used as the initial source of class or 
subclass statistics . 

B. Test and/or training fields will be used to evaluate classi- 
fier performance . 

C. DO fields will automatically be identified as the categoiy 
"other," regardless of pixel data values, except when process- 
ing test fields, or training fields for the category "Other." 
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4. 1.1. 3.1 
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D. Pixels lying in DU fields will automatically "be identified as 
"belonging to the category Unidentifiable, regardless of the 
pixel data values, with the following two exceptions. Pixels 
lying in both DO and DU fields will be identified as belonging 
to the category "Other." When classifying test fields, or 
training fields for the category "Unidentifiable", all pixels 
will be classified on the basis of pixel data values . 

E. MLA fields will be used in an algorithm (below) to compensate 
for atmospheric variations. 

4 . 1.1. 3. 2 Field definition: The system shall generate a field verification 

overlay at the same scale as LACIE AI imagery. 

Field definitions to be used in further processing shall be obtain- 
able by each of the following routes : 

A. 1. By input of coordinates selected externally to the system. 

2. By output and reinput of coordinates selected interactively 

on the system. 

B. Interactively by online field definition. In this mode, the 
analyst may request recall of the stored coordinates as a 
starting point. The interactive operation, however, will not 
update the fields data base. 

Fields can contain any number of pixels, pixels/line, or lines per 
column from one up to the limit set by segment dimensions • Fields 
may have from 2 to 10 vertices. Minimum field width shall be one 
pixel . 

it^.1.1.3.3 A ’replace’ option should exist for the Fields Data Base. The 

’replace' option would allow existing fields and statistics to be 
deleted from one category, and in a single operation, then be in- 
cluded in a different category. 

i^.l.l.3.i^ The system shall generate an overlay of current dot labels at the 
same scale as LACIE AI Imagery. Tbe overlays shall be on the same 
scale as the PFC Product 1. The overlay shall be transparent with 
each dot designated by a unique color or gray shade by category. 

This color or gray shade shall be consistent between segments . The 
overlay should carry the data of the update and the segment identi- 
fication. 

Interactively, dots, labels and types to be used in further process- 
ing shall be obtainable by each of the following routes: 

A. By input of coordinates, labels and types selected externally 
to the system. 

B. By online dot definition. Analyst may request recall of stored 
dots as starting point. 

C. Dots with labels and types shall be displayed as an overlay to 
the television monitor, either color or gray level. Present 
labels are W, W, X, DO, and DU. Dot types are S for starting 
vector, L for labeling vector or B for bias correction vector. 



D. Dot definition changes shall be controlled by analyst via key- 
board inputs. Modification to the dot data base by analyst 
shoiild be restricted by use of read/write keys to insure data 
base security. 

I 

4. 1.1. 4 Computation of statistics .- The system shall be capable of aggre- 
gating the training fields from up to three segments for all further 
processing defined in this section. That is, the system shall be 
able to process these fields as though they were located in a single 
segment. The system shall also be capable of merging statistics 
obtained in up to three nonagregated segments. 

Three modes shall be available for computing statistics from train- 
ing data as follows : 

Mode 1 : Input training fields are used to directly define class ’ or 

subclass statistics without further operations. This mode 
is the defavilt option. 

Mode 2a: Input training fields are used to provide starting 

vectors and data points for clustering algorithm. Clus- 
• taring algorithm develops clusters. Cluster chaining 
algorithm is used to chain similar clusters together. 

The (chained) clusters are then used to compute subclass 
statistics . 

Mode 2b; As above without chaining of clusters. 

Hard copy output products from the clustering modes will consist 
of the following: 

A. Cluster Map: Two types of color coded cluster maps shall be 

produced in transparency form at the same scale as other seg- 
ment imagery (ref. sec. 5»1.6.l) whenever a segment is batch 
processed. These cluster maps shall provide the capability to 
display all the clusters produced in any single clustering run. 

A.l An \mconditional cluster map 'shall be produced in which 
the pixels belonging to each cluster are displayed as a 
different color where the colors are chosen in order from 
a color table . A capability shall be provided to input 
a default color table. 

A.2 A conditional cluster map shall be produced in which all 
unconditional clusters with the same categorj^ label are 
displayed as the same color but with a different color for 
each category. In addition, each conditional cluster shall 
be displayed using a different color. A conditional clus- 
ter is a cluster whose LI distance from the subclass or 
-dot used to label it to the cluster mean is- greater than a 
user input threshold. The colors used to produce this map 
shall be selected in order from a color table. A capabil- 
ity will exist to input a default color table. 



In order to^ aid the analyst in reviewing the cluster maps produced 
by a hatch job, each map shall be annotated with a color key in the 
form of an example of each color that was used to make the map in 
the order in which the colors appear in the color table* The 
color key for a cluster map shall appear adjacent to the map i*e* 
between the image and the header (ref. sec. 5*l-6*2.6 - Image 
Header) . 

Cluster maps shall indicate the boundaries and identities of all 
fields by a technique which does not prevent the user from ascer- 
taining the cluster of any pixels. It is permissible to use mxiLti- 
ple sheets to accomplish this. 

B. Reports comparable to the ERIPS detailed clustering reports 
including options used in clustering, cluster-class matches, 
cluster population and cluster names. 

C. Distance tables including intra- and interclass distances, for 
all classes included in the clustering run. 

D. A report for each cluster generated which shows categories 
assigned to that cluster, the labeling dot closest to that 
cluster, and a listing of individual dots assigned to the 
cluster with the distance from each dot to the cluster mean. 

The report shall contain the information in the following 
format : 

Cluster Name ; Labeling Dot ; Categories ; 

DOTS: DISTANCE: DOTS: DISTANCE : 

XX, X XX. X 



XX. X XX. X 

Mode 3: Input training fields are assumed to represent clusters. 

The Bhattacharyya distance measure is used to decide 
whether or not to chain fields together. The (chained) 
fields are then used to compute subclass statistics. 

Hard copy output products from Mode 3 shall consist of the follow- 
ing (menu to be furnished): 

A. Intraclass distance tables H[(i,j), (i,k)] for all fxelds in 
each class (see flow charts for details). 

B. List of chained fields. Number of pixels per chain. 

The modes shall be individually selectable by class by segment. 
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The clustering algorithm shall he an interactive algorithm eq,uiva- 
lent to ISOCLS, and shall process up to l6 channels of data, or any 
subset of channels defined by the tiser. The user shall be able to 
define such subsets on a per segment, or per aggregation of seg- 
ments basis. All segments to be classified as an aggregation will 
have the same channel subset. 

A. The iterative and one-pass clustering algorithms shall each 
provide a capability to generate up to 60 clusters whose sta- 
tistics can be saved for subsequent use in classification. 

B. The iterative clustering algorithm shall have the option of 
acquiring starting vectors from the one-pass algorithm in batch 
as well as interactive modes while retaining the current op- 
tions. Provision for batch control of the parameters of the 
one-pass algorithm shall be made. 

C. The option shall exist to use either the L-1 or the L-2 
distance for assignment of pixels to clusters. Specifi- 
cally, these distances are: 

1. (L-1) D1(X,Y) = ^ lx - y I , and 

p=l ' P p' 

2. (L-2) D2(X,Y) = T (x - X )^ 

- ■ ' p=l p p 

* - ' 

' where n is the number of channels, x is the vector 
{x^, x^, x^}, and y is the vector 

{y^, 72» •••» 

D. On a split iteration, the iterative clustering algorithm shall 
split clusters in order of their maximum standard deviation 
from largest to smallest. On a combine iteration, the itera- 
tive clustering algorithm shall“Combine pairs of clusters in 
order of the distance between them from smallest to largest. 

Any cluster can be split or combined only once per iteration. 

E. The option shall exist in clustering to correct the dynamic 
range of each Landsat acquisition in a composed image for the 
effects of sun angle, with either the L-1 or the L-2 distance. 
This correction is done by altering the distances with a 
diagnoal matrix A. The new distances are then 

DIS (x,y) = Dl(Ax,Ay) = A^^ |x^ - y^j 

and 

D2S (x,y) = D2(Ax,Ay) - A^^ (x^ - y^)^ 

P. Upon completion of the last iteration of the iterative cluster- 
ing algorithm, all clusters with less than P(K) = n + 1 (where 
n - i channels) pixels shall be eliminated and those pixels 
assigned to the remaining clusters . 

G. The number of sets of default values for clustering parameters 
shall be increased from two to eight or more (with analyst op- 
tion to override for each parameter) . The strip generation 
parameter and the select statistics box parameter shall be 
removed from and the E2 parameter shall be added to the list of 
default parameters. 
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H. The option to eliminate DO (designated other) and/or DU 
(designated unidentifiahle') fields from clustering is required. 
The option to input maximum and minimum radiance values in each 
channel and to exclude from clustering those pixels outside the 
interval (min, max) is also required. 

I. There shall he an option to label each of the final clusters 
with the category of the closest starting or labeling vector. 
The distance used shall be the same as that used in pixel 
assignment. There shall be an option in batch processing 

to input a label for a given cluster number or to delete a 
specific cluster's statistics from classification. 

J. There shall be an option to input a cluster distance threshold 

for use in Subclass Match. If any cluster is further from the 
nearest subclass than this threshold distance, a conditional 
cluster map will be generated in which each such cluster is 
assigned a color that indicates: (l) which subclass it was 

matched to, (2) that it exceeded the threshold distance, and 
(3) that a conditional cluster map will be generated every 
time a segnent is clustered. 

K. The cluster chaining algorithm indicated in the flow charts 
shall be used to chain similar clusters together. 

L. The user shall have the capability to select arbitrary subsets 
of training segment channels for further processing, or to 
allow an automatic feature selection algorithm to select an 
optimum channel subset. Furthermore, he shall have the capa- 
bility to associate any acquisition date of training segment 
data with any biological phase for classification, and he shall 
have the capability to utilize training data from any one of 
three segments (s) to classify any other segment. 

M. The user shall have the option to compute bias corrected cate- 
gory percentages from-cluster populations. This correction 
shall be made automatically in the event that all clusters are 
given the same category label. Correction involves. 


p = the number of categories to be bias corrected. 

K = the number of pixels excluded from clustering by the min/ 
max, radiance test. 

= the number of pixels in DO fields 

a^ = number of pixels assigned to cluster number i, 

s = {set of clusters that were labeled category C } 
q ° q 

N = E a. = Nuniber of pixels assigned to all clusters 
^ ies ^ labeled C , 

q. q. 

s = {set of all clusters with labels other than A} 

IT = E a. = Humber of pixels assigned to all clusters not 
ies labeled A. 

B = N + HS 

= the number of pixels assigned to the totality of clusters 


m 


q.>r 


Ug 5 . « . C 

= The nimiber of bias correction dots labeled C 


whose label is other than A or C^, >^2 5 p 


but 


assigned to the totality of clusters labeled C for q and 
r < p. 


U-T 



m ., = The number of bias correction dots labeled C but 
1>P+1 q. 

assigned to the totality of clusters whose label is 
other than A or for t ^ p. 


m 


P 

S 

< 1=1 


m 




+ m 


p+l,r 


= The number of bias correction dots 

assigned to the totality of clusters 

labeled C . 

r 


m 


P+1 


The nmber of bias correction dots 
totality of clusters labeled other 
q. ^ P- 


assigned to the 
than A or C for 

q 


m 

a = —3^ or a = 0 if m - 0 for q. ^ p and r ^ p+1 
q,r m^ g,r r i' j- 


P 




100 


For each 
is 


q ^ Pj P+l-limensional vector whose 


th 


component 


Y 


'k l2 (a \ (l-a 

= ^ xiOO i 

q^r B m -1 

L J r 


-o ') 


VAR 




for p 


N 

B = X 100 -PC g - P 
q B q ^ ^ 

The following segment summary information shall be output from 
cluster bias correction to the dot summary report. 

P^C^^ through 

1 p’ p+1 

^ for g < p, r £ p+1 

g>x 

V for q ^ p, r < p+1 

qjT 

B for q ^ p 

q 
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VAR 


|^P(C^)J for q i p 


and 

\ • V f ^ ^ • V ['’('o] 

where Z = 1.6^5 

U.1.1.5 Modification of statistics .- The capability to perform the follow- 
ing operations on class or subclass statistics shall exist: 

A. Sun an^e corrections 

B. Mean level corrections 

The user shall he able to select from the following options : 

A. No corrections. (This shall be the default option.) 

B. Sun angle correction only. 

C. Mean level adjustment only. 

D. Sun angle correction followed by mean level adjustment. 

These options shall be available to the user on a per segment to 
be classified basis . Class and subclass statistics shall be 
available as a user optional output following the optional modifica- 
tion of statistics. In classifying a merged image, the user shall 
have the option to classify the training segments using unaltered 
statistics . 

The output shall be available in the following forms : 

A. Line printer, microfiche, or microfilm. This output shall in- 
clude means, covariance“matrices , and correlation matrices. 

B. Magnetic tape. It shall be possible to input statistics at 
this point by use of magnetic tape. The above capabilities are 
operational. 

4. 1.1. 5.1 Sun angle correction: These corrections shall be performed using 

coefficients obtained from a table organized by sun angle of acqui- 
sition of training segment versus sun angle of the equivalent acqui- 
sition of the recognition segment. Note that, since up to four 
different acquisitions of each segment are combined to create up to 
l6 channels of data prior to beginning pattern recognition, one to 
four pairs (of sion angle of training segment data acquisition versus 
sim angle of recognition segment data acquisition) exist. For each 
such pair a multiplicative and an additive correction exists for 
each Landsat band. 

These corrections are represented by 4 x 4 diagonal matrices. A, 
and 4-element vectors , B . An A and a B exist for each sun 
angle of acquisition pair (i.e., SEA FT and SEA PR), and are denoted 

A (SEA FT, SEA PR), B (SEA FT, SEA PR) 
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where 


SEA PT = Solar elevation angle of acquisition P (P=l-U) of training 
segment 

SEA PR = Solar elevation angle of acquisition P (P=l-U) of recogni- 
txon segment 

All off-diagonal elements of A are zero. The on-diagonal elements, 
a. . are decimal numbers X.XXX between 0 and 5 ♦000. 

XI 

The elements of B. 

X 

b^ are decimal numbers XX. XX lying between — 99«99 and 99* 99* 

i is a subscript indicating the Landsat spectral band for which 
the a..’s andb.'s are valid. 

XX X 

The values of A and B shall be stored in a table extending 
from SEA = 20® to 70° for both training and recognition segment. 

The table shall be addressable at increments of 1®, although it is 
not necessary to store values of A and B in the table incre- 
ments closer than 5°- Interpolation may be used within a 5° 
interval . 


For a given data set, consisting of acquisitions M, N, 0, P, 
form a sun angle correction matrix 



/a (sea MT, 

SEA 

MR) 

(0) 

(0) 

(’0) \ 


■ (0) 

A 

(SEA 

NT, 

SEA NR) 

(0) 

(0) 1 

AT = 









1 (0) 

(0) 



A (SEA or, SEA OR) 

1 


\(o) 

(0) 



(0) 

A (SEA PT 

, SEA PR) / 

where 

each 

entry 

is a 

4 X 

4 matrix. 

Also form 

a sun angl( 


vector 


BT = (B (SEA MT, SEA MR) , B (SEA HT, SEA NR), 
B(SEA or, SEA OR), B (SEA PT, SEA PR)) 
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Then for each^ suhclass, the sun angle corrected mean vector, yc, 
is formed from the original mean vector y as follows : 

yc = y X AT + BT 

The s-un angle corrected covariance matrix Fc is formed from the 
original covariance matrix F as follows: 


Fc = AT X F X AT. 


The above capability is operational. 


1^.1. 1 . 5 *2 Mean level adjustment: The system shall compute. the mean vector in 

the training segments and in the segment to be classified over the 
MLA field(s) in each. If sun angle correction was requested, the 
system shall correct the MLA vector for the training segment as in 

the previous section. Then compute^ Ay = y (training segment) - y 
(segment to be classified). Then for each class or subclass mean, 

y (class), compute the mean level corrected vector y (class) as 

follows: y (class) = y (class) - Ay. The above fiinctions are 

operational. 

14.1.1.6 Feature selection requirements . - The system must provide an optional 
feature reduction procedure employing up to and including 16 channels 
(The defatilt option is to by-pass feature reduction. ) The separability 
measure criterion shall be the Bhattacharyya distance measiire. The 
featvire reduction algorithm must be able to process statistics for all 
the classes and subclasses defined within the data base for each segment. 


The featTire reduction procedure must be able to operate in the 
following modes; 


Mode 1: This mode is the default option if the number of channels 

is over ei^t . Optimal channel subsets will be selected 
by a "without replacement" procedure xxntil the subset sep- 
arability is at least x percent of the total channel 
separability. As a user option, black-and-white (B&W) 
film may be produced for best single channel and/or color 
film may be produced for the best three channels . 

Mode 2 : This mode uses all the channels and is the default option 

if the number of channels is eight or less . 


Mode 3: This mode uses a subset selected by the analyst. 
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Mode 4: This mode will compute Bhattacharyya disteince separability 

for acquisitions and combinations of acquisitions. 

A priori probabilities for use in the computation of Bhattacharyya 
distance shall be ret\irned from the data base \anless overriden by 
analyst input. Override may be accomplished in both batch and 
interactive modes . 

Hard copy output products shall consist of the following: 

A. Number and identification of channels and the total number of 
channels available. 

B. Interclass, intraclass, other, and total separability as -defined 
in the flow charts computed using all the channels and the 
extracted set of features. 

C. The ratios of 1 of the above computed using all the channels 
and the extracted or analyst selected set of features. 

4, 1.1. 7 Classification .- The classifier will accept or retrieve the follow- 
ing inputs for each segment to be classified. This capability is 
operational. 

A. Class or subclass statistics for each class (as modified by the 
feature extraction process). 

^ priori probabilities for each category for the segment to be 
classified will be retrieved from the fields data base. The 
user will have the capability to override this value through 
category input or by selecting computation of proportional siib- 
class ^ prioris . 

C. Threshold values for each category for the segment to be classi- 
fied will be retrieved from the fields data base. 

D. DO and DU field definitions will be r etrieved from the fields 
data base. 

E. Test and training field definitions will be retrieved from the 
fields data base. 

F. Channels to be used. 

G. Pixels from the segment to be classified. 

H. Thresholding should be implemented on a subclass basis with 
analyst-selected percent thresholding down to 0.1 percent in 
increments of 0.1 percent. Based on the number of pixels in 
the subclass training field and the nianber of channels tised for 
classification, the computer will automatically determine the 
required numerical input assuming a Fisher's F distribution. 

The default threshold value should be selected by the computer 
and be based on the single numerical entry per category which 
gives the most acc\irate threshold next to the single nvanerical 
entry per category which gives the most accurate threshold next 
to the single numerical entry per subclass model, i.e., a de- 
fault value per category level which does not significantly 
change the level of thresholding. 

The analyst should have the option to bypass the automatic 
thresholding and input values at the subclass level. 
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The threshold' option coiild he used at any of the following 
entry points: “ (l) Class summary menu;' (2) 'all classification 
image menus; (3) fields data hase card interpreter; (i^) process 
control card input . 

Ij.l.l.T.l Two types of classification will he performed as follows: 

A. Classification of training and test fields. 

B. Classification of segments for acreage estimation. 

In the classification of test and training fields, pixels falling 
in DO and DU fields will he classified as though they were ordinary 
pixels. In classification for acreage estimation, a.i 1 points in DO - 
fields will he classified as other and all points in DU fields will 
he classified as imidentifiahle , unless they also 'lie in DO fields. 
Classification will assign pixel to a certain subclass and record 
that assignment in the Classification Work Data Set' (Appendix B-13) 
along with an associated likelihood value. The density function to 
he used in the assignment decision involves the following contribu- 
tion from a given subclass : 


P (x,k,i,j) 




- 1/2 


e 




where x 

k 

i 

j 

a 

"k 

m. 

1 

d 


is “vector representation of the pixel to be 
classified 

is the category identifier 
is the class identifier 
is the subclass identifier 

is the priori probability corresponding to cate- 

gory k 

is the number of classes in category k 

is the number of subclasses in class i of cate- 

gory k 

is the dimension of pixel vector x 




k ,i , j 


is the mean vector of subclass (k,i,j) 
is the covariance matrix of subclass (k,i,j) 
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U. 1.1. 7.2 


In th.e nominal default classification to the category level, the 
following sum would be computed for each categoiy k: 


“k 


m. 

1 


P(x,k) = 



P(x,k,i,j) 


i=l j=l 


The pixel ¥Oiald then he assigned to the category with the largest 
associated sum. Specifically, it shall he assigned to the subclass 
within that category which has the maximal P(x,k,i,3). Stored with 
that subclass assignment shall he the following associated likeli- 
hood value to he used for thresholding: 



If this number is larger than the chi-sguare value, the pixel is 
classified into the category threshold. 

For any category k for which classification to the class level 
has been specified^ sums of the form 


m. 

1 

P (x,k,i) = ^ P (x,k,i,j) 

j=l 


would enter the overall consideration for largest density sum as 
above. The specific subclass assignment and associated likelihood 
value would be derived in the same manner. 

The classifier will have the following outputs for each recognition 
segment : 

A. Classification Map 

1. The classification map shall indicate wheat versus other 
versus thresholded versus unidentifiable, in a transparency 
form, at the same scale as LACIE PI imagery for the recogni- 
tion segment alone. 

2. The classification map shall indicate the boundaries and 
identities of all fields by a technique which does not pre- 
vent the user from ascertaining the classification of any 
pixels. It is permissible to use multiple sheets to accom- 
plish this. 


It-l4 



B. 1. 


1 . 

2 . 


D. 


E. 


F. 


The grey shades or colors used in the classification maps 
shall he consistent by category between se^ents , yet re- 
taining contrast. For example, wheat on Segment A should 
be represented the same as wheat on Segment B. CAMS shall 
have the capability to specify default values. 

Threshold values used and actual number of pixels 
thresholded/class , category and subclass. 

Date and user identification of input of threshold values 
used. 

A priori probabilities used per category. 

Date and user identification of input of ^ priori values 
used. 

Classes and subclasses used, source of each subclass, i.e., 
clustering, B-chaining, etc. 

(Transformed) band means and variances for each class or 
subclass . 


Classified percentages for each category C^, C^,. 
bias corrected where : 


C to be 
P 


C^ = name of category to be bias corrected, 1 5 q ^ p. 
p = number of categories to be bias corrected. 


= mmiber of pixels in DU fields. 

= number of pixels thresholded in classification 
= number of pixels classified as Categoiy A 

B = 22,932 - \ 

- number of pixels in all DO fields 

A = the name of a category vhose pixels are to he excluded 
from classification percentages. A shall he a user input 
"with default A=X. 

X - usually an unidentifiable category (clouds /cloud shadows). 


= the number of pixels classified in category C^. 


N 


p+1 


= the mnnber of pixels classified into a category other 


than the categories A, C- , C, 


N 

= ^ X 100 ; 


^p+1 


- 

B 


1 ’ '" 2 ’' 
X 100; 


. C, 


W 

^DO = ^ X 100 . 

JJ 


n 1 Number of pixels in DO fields 

196*117 ^ number of pixels in DU 
fields - number of pixels in category 
unidentified 
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I 


2. Corrected percentages for each category 
For an index, r, 1 r ^ p+1 


P C 


p+1 N 


a I X 100 

q,r B 


r=l 


where 


m 




number of bias correction dots labeled C by the 
analyst and by the classifier 


m = number of bias correction dots labeled C by the 

g,p+l q 

analyst and as some category other than A, 

, . . . , C by the classifier . 

2’ ’ p 


m 

r 


Nimiber of bias correction dots labeled C by the 
classifier. ^ 


m = number of bias correction dots that were labeled by 
^ the classifier as some category other than A, , 

• » P* 

The summary shall be output such that a sequential listing 
of type 1 dots is output first, followed by all dots (other 
than type l) classified C^. 

Note, any bias correction dots that were thresholded in 
classification or that fall in DU areas shall be excluded 
from m 


a 


m 

<l>r 



0 


\,p+l/“p+l 


if one or more bias correction dots 
were classified C for q=l,2..., p and 
r=l,2, . . . , p 

if no bias correction dots were classi- 
fied C i.e. m =0 
r r 

for q=l,2 , . . . , p and r=p+l 


A dot processing summary shall be provided for each segment 
processed (whether classification is performed or not). 

In general, this summary shall contain the following infor- 
mation sequentially listed with respect to dot numbers : 

Dot number 

Dot pixel, row ntimber 
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Dot classification category {as detjermined by classifier if 
classified) 

Dot cluster label (label of cluster to -which dot was 
assigned by clustering algorithm) 

Dot category label (if specified by analyst) 

Dot type (S-starting vector, L-labeling vector, B-bias cor- 
rection, DO, DU) 

Dot green number 
Dot bri^tness 

The summary shall be output such that a sequential listing 
of starting and labeling dots is output first, followed by 
all dots (other than starting or labeling dots) classified 
then Cg sequentially through C^, followed by dots 

classified other than. A, C^, Cp,..., C , followed by DU, 
threshold, then DO. 

In addition, the following segment summary information shall 
be output to the dot s-ummary report . 

P(C^) through P(C ) 

I P 

Np, 

through ^^p+1’ 

-through 

“ 2.1 “2,ptl 

through 


a > through a , - , 

p,i ® p,p+i 



N 

B = :^ X 100 - P/C ) 

q B \q.) 


U-IT 



and 


- Z ^VAE + zjvjm 

where z ^ I.6U5. ' 

H. ^Timber of -pixels in DU fields 

196*117 

I. Source of training data: 

1. Training segments(s), date(s) of acquisition. 

2 . Training fields . 

3. External statistics. 

J . Channels used . 

K. Channel map used to relate training segment channels to recog- 
nition segment channels. 

L. Identity of recognition segnent. 

M. Acquisition dates of recognition and training segments, includ- 
ing all flaggings of marginal data. 

N. Whether nin was semiautomatic or interactive. 

O. Playback of analyst semiautomatic mode inputs. 

P. Playback of analyst interactive mode inputs. 

Q. Identifying run number and date of run. This item shall appear 
on all individual products. 

R . Biological phase . 

S. A classification summary similar to the ERIPS classification 
summary, generalized to include subclasses, pixels classified 
into each subclass, pixels threshold for each subclass, for all 
test and training fields in recognition and training segment(s). 

T. As item S, except by class, and category. 

The following reports shall be included with the classification 

report when applicable: 

U. A fields definition report similar to that generated by the 
ERIPS, but including subclass name for each field. 

V. A clustering report including channels used for clustering, 
cluster population, class matches, cluster names, a detailed 
clustering report , interdistance clustering report , intradis- 
tance clustering report, cluster parameters, definition of 
cluster maps generated for PEG and also output on microfiche, 

W. A chaining report including fields /clusters chained, an intra- 
distance report (B-distance; initial and final), initial number 
of subclasses and subclass identities, final number of subclass- 
es and subclass identities, number of pixels per chain. 

X. A feature selection report including channels available, chan- 
nels selected, interclass separability, total separability with 
all channels , total separability with selected subset of 
channels . 

y. User name • 

Z. For Interactive runs, if interactive field changes are made, a 
complete definition of fields , ^ priori probabilities , and 
threshold values used for classification shall be reported. A 
field boundary overlay shall also be generated and referenced. 
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AA. Diagnostics and incorrect inputs-. 

BB. The capability should exist to print selected pattern recogni- 
tion reports. 

4. 1.1. 7 *3 Classification output products shall he distributed as follows: 

A. One copy, located on microfishe, of each of items A. 2 to BB 
above will be delivered to EOD for the data processing analyst 
associated with the recognition segment. 

B. Items A.l, B.l, C.l, F, G, H, I.l, 1.3, K, L, M, N, Q, R, S, T, 
U, V, W, X, Y, Z, AA, and BB shall be stored in the ISERS for 
each classification run; later runs shall not supersede earlier 
runs . 

C. The following parameters shall be stored on magnetic tape for 
each (batch) classification run. 

For the recognition segment: 

1. Sample segment number 

2 . DPR number 

3. RTxn ID 

4. Bhattacharyya separability 

a. All channels 

b. Channel subset used 

5 . External statistics tape number 

6. Rumber of recognition segment subclasses 
T. Number of merged image fields 

8. Number of pixels in DO fields 

9. Number of pixels in DU fields 

10. Date of classification 

11. Number of clustering r\ms 

12. Number of chaining runs 

13. Nimiber of acquisitions used for classification 

14. Channels used in classification 

15. Nominal biological windows 

16. Segment type identifier ( spring/winter ) 

17. Sample segment use indicator (training or not training) 

18. Dates of acquisition on merged image 

19. Segment number, segment type indicator, and acquisition 
dates for each training segment 

20. Threshold and a priori values used in classification (by 
category) 

21. Subclass clissification resiolts 

a. Subclass name 

b. Statistics manipulation flags 

c. Pixels classified into subclass 

d. Pixels threshold from subclass 

For each field in segment : 

22 . Field name 

23. Subclass name 

24 . Field type 

25 . Field population 

26 . Number of subclasses 

27. Name of each subclass 
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28. Mean level adjustment. Sun angle correction flags for 
that (sul) class 

29. Pixels classified into that subclass 

30. Pixels threshold from that subclass 

31. Means and standard deviations for each class 

32. Information from items U. 1.1. 7. 2 F, G, H, and Q 

33. Cluster match and distance information 

This tape shall be stored by ISERS in the physical library for 
use by CAMS and CAS. 

3U. Acquisition separabilities 

35. Cluster option indicators 

36. Cluster/dot report 

37. Cluster bias correction information 

38. Cluster green number & bri^tness value for each acquisi- 
tion used in clustering. 

The analyst evaluation of the recognition segment shall be manually 
passed to CAS and stored in the physical, library by ISRRS. The 
following evaluation codes shall be used: 

Code 01: Not machine processed - clouds, haze, snow, etc. 

Means that the segment cannot be processed through the 
system because clouds, haze, etc., make interpretation 
and analysis impossible. 

Code 02: Not machine processed - confusion crops or other inter- 

pretation problems . 

This code should be used when a segment cannot be proc- 
essed because of interpretation difficulties , especially 
when confusion between wheat/small grains and other crops 
is such that a wheat /smadl grains estimate cannot be de- 
termined. The rules for such a decision are to be nego- 
tiated between CAMS and CAS and included in the CAMS De- 
tailed Procedures . 

Code 03: Bad data - due to technical problems, not reordered. 

In this case the segment cannot be processed due to 
technical problems arising from unsatisfactory histo- 
gram, excessive scanline drop, etc. If the segment is 
not reordered, this code is used. 

NOTE: If the segment is reordered, no code should 

be listed and no data sheet should be passed 
to CAS since this segment will again be sent 
through the CAMS system for evaluation. Any 
segment that cannot be processed due to clouds , 
etc., or technical problems, should be passed 
as 01 or 03. 
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Code 05: Wot machine processed - dormancy. 

In cases where recognition is a problem due to the 
crop being in a state of dormancy this code is used. 

Code OJ: Wot machine processed - preemergence. 

In cases where the acq,xiisition is prior to the criteria 
established for fall wheat recognition (TBD) this code 
is used. 

Code 09: Wot machine processed - multiple acquisitions. 

Code 09 was developed to take care of the problem of 
processing multiple acquisitions of a segment at the 
same time. When multiple acquisitions are available at 
one time , only one segment , which is detennined to be 
the better acquisition for interpretation and analysis 
by the CAMS analyst, is selected for complete processing. 
Other acquisitions are listed as Code 09 to indicate that 
they had been reviewed but were not processed through the 
system. This code enables CAS to review the segment re- 
sults and account for each acquisition and how it was 
evaluated by CAMS. 

Code 10: Unsatisfactory - unsatisfactory results for segment. 

This code is to be used for any acquisition that proc- 
essed through the system and, based on CAMS evaluation 
procedures, designated unsatisfactory. 

Code 12: Unsatisfactory - no significant change. 

This code is used when the new acquisition is evaluated 
to have no significant change from the previous unsatis- 
factory evaluation for the segment. 

Code lU: Unsatisfactory - rework, reevaluated segment. 

This code is used when a segment which was previously 
passed to CAS is reworked. 

Code 18 : Unsatisfactory - machine processed multitemporal 

analysis . 

This code is used when more than one acquisition date is 
used to produce an unsatisfactory proportion estimate. 

WOTE: All acquisition dates used in processing should be 

listed on the CAl'IS Evaluation Form. 
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Code 


Code 


Code 


Code 


Code 


Code 


Code 


Code 


20: Marginal - marginal results for segment. 

This code is to be used for any acquisition that is 
processed through the system and, based on CAMS evalua- 
tion procedures, designated marginal., 

22; Marginal - no significant change. 

This code is to be used when the new acquisition is 
evaluated to have no significant change from the previous 
marginal evaluation for the segment. 

2k: Marginal - rework, reevaluated segment. 

This code is used when a segment which was previously 
passed to CAS is reworked. 

28; Marginal - machine processed multitemporal analysis. 

This code is used when more than one acquisition date is 
used to produce a marginal proportion estimate, 

NOTE: All acquisition dates used in processing should 

be listed on the CAMS Evaluation Form. 

30: Satisfactory - satisfactory results for segment. 

Ihis code is to be used for any acquisition that is 
processed through the system and, based on CAMS evalua- 
tion procedures, designated satisfactory. 


32: Satisfactory - no significant change. 

This code is to be used when the new acquisition is 
evaluated to have no significant change from the pre- 
vious satisfactory evaluation for the segment. 

34: Satisfactory - rework, reevaluated segment. 

This code is used when a segment which was previously 
passed to CAS is reworked. 

36: Satisfactory - less than manually (hand) counted. 

This code is to be used for any segment in which the pro- 
portion estimate is manually counted rather than machine 
processed. 

NOTE: This category would not be used if the segment was 

a rework segment. Code 34 should be used. 

Segments in an agricultural area that have a 0% 
proportion estimate are to be designated Code 36. 


4-22 



Code 38: Satisfactory - machine processed multitemporal analysis. 


This code is used when more than one acquisition date is 
used to produce a satisfactory proportion estimate. 

Code UO: Segment is totally nonagricultural . 

This code is used when the segment is evaluated as hav- 
ing no agriculture at all, i.e., no discernahle field 
patterns . 

NOTE: Se^ents in an agricultural area that have a 0 % 

proportion estimate are to he designated Code 36. 

The ISRRS shall he designed to store a total of U0,000 such reports; 
the storage available for any given segment shall he limited only hy 
total available storage space. 

Storage of TBD portions of TBD number of these reports for TBD 
period following end of crop year is required- (Category 3) . 

The above results shall be accessible by recognition segment number 
and by training segment number (s) . 

D. A paper printout of a report, in the format as described in 

Appendix A, will be delivered to EOD for the data processing 

analyst associated with the recognition segment. 

4. 1.1. 8 gpjTni automatic mode .- The referenced flow charts indxcate the semi- 
automatic mode. The semiautomatic mode shall, as a minimum, support 
the following operations: 

A. Preprocessing.- An option to produce spectral plots at C&I 

loading (Ref. Sec. U.1.5). 

B. Classiciation.- The inputs shall be the following: 

1. A collection of one or more recognition segments for class- 
fication processing - no default. 

2. A set of from one to three training segments, or external 
statistics for the entire collection of recognition seg- 
ments defined in (l) . 

3. Channel map(s) relating channels of each acquisition of 
designated training segment (s) to channels 1 through 16. 

4. Channel subset instructions for processing entire collec- 
tion of recognition segments . This channel -subset must 
exist for each recognition segment and training segment 
defined in (l) and (2) above. 

5. Statistics computation mode by class and by training 
segment . 

6. Instruction to perform sun angle correction for selected 
recognition segments . Default- shall be no correction . 

7 . Instruction to perform mean level adjustment for selected 
recognition segments. Default shall be no adjustment. 

8. Instructions to classify test and training fields in 
training segments with unaltered statistics. 
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9. Featttre reduction laode(s). The default is to hy-pass 
feature selection. 

10. User name. 

C. Computation of statistics for external output. Inputs are 
identical output. Inputs are identical to (A2) through (AlO). 

D. Updates of field definitions, a priori prohahilities . thresh- 
hold values. 

E. An. option to produce spectral -plots (ref. SEC. 4.1.5). 

t 

The above capabilities are optional. 

4. 1.1. 9 Interactive mode .- The interactive mode is viewed as a facility 
for processing ‘those segments for which satisfactory results can- 
not "be achieved in the semiautomatic- mode . It shall combine the 
features of the semiautomatic mode and the current ERIPS. 

Override of previous results may be accomplished by analyst addi- 
tion or redefinition of fields, changes in classification, preproc- 
essing modifications to statistics, selection of channels, etc. 

This shall be accomplished by online interaction. 

Interfaces between the interactive and semiautomatic modes shall 
be as follows : 

A. - Access to a common data base. (Some means will be .required to 

minimize impact of the semiautomatic mode to the interactive 
mode and -vice "Tersa when the two modes require access to two 
segments which may not be simultaneously accessible.) 

B. Statistics output from the interactive mode may be used as an 
input to the semiautomatic classifier. 

C. Output from the interactive mode may be reinput to update the 
ISRRS. The ISRRS shall indicate the mode of update in its , 
records . 

The interactive mode is operational. 

4.1.1.10 Selection of training segaents .- The system shall have a mode by' 
which it assists the analyst in specifying the best training segment. 

For manual training segment selection, access to the following is 
required; 

A. Acquisition status at JSC of all segments. 

B. Acquisition status at GSPQ of all segments. 

C. Nominal acquisition dates for each segment. 

4.1.2 Training 

(Req-’d by: CAMS; Cat. 1; Impl. Resp. : CAMS/LACIE AVSB) ‘ 

It is required that a training program adequate to indoctrinate the 
available personnel with a broad understanding of remote sensing, 
wheat growing, and recognition processing, together with a detailed 
knowledge of the LACIE CAMS data system and the procedures for its 
use be planned and operated. 
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4.1.3 


4. 1.3.1 


Procedures Development 

(Reg'd by CAMS: Cat. 1;' Impl. R'esp. 


CAMS/LACIE ASYB')- 


4. 1.3. 2 


4. 1.3. 3 


4. 1.3. 4 

4. 1.3:5 
4. 1.3. 6 

4.1.4 


Development of the following procedures is required. 

> ^ I ' 

Image interpretation procedure A procedure exists which allows a 
trained analyst interpreter, using imagery produced from Landsat 
MSS data, to: 

A. Determine and categorize -those ground cover classes which can 
be distinguished from wheat and those which cannot. 

B. Identify training fields for each of those classes by line 
and pixel coordinates . 

C. Label those training fields in a manner suitable for entry into 
the ADP system. 

D. Perform the necessary quality control associated with the above 
operations. 

Signature extension partitioning procedure .- A procedure exists 
which allows a trained data analyst, using information characteriz- 
ing local soil types, meteorology, farming practices, crop calendars 
and Landsat MSS data, to divide wheat growing regions into regions 
of nearly constant wheat signature and to associate regions of 
similar signatxn’es , to designate a set of training segments , and 
to redesignate training segments as required. 

Training segment selection procedure .- A proced^^re exists to govern 
the selection of a training segment to be utilized in the processing 
of a recognition segment, given the above signature extension parti- 
tioning, the acquisition histories of the involved segments, and 
recent meteorological data. 


Automatic data processing procedure .- A procedure exists to govern 
the choice of options used in the classification of segments given 
the recognition segment, training segment, acquisition histories 
and ancillary data. 

ADP evaluation procedtire .- A procedure exists to be used in evalu- 
ating the results of ADP classifications to assess their accuracy. 

Rework procedure .- A procediire exists to be used in reworking 
classification results of inadequate accuracy. 

Green Umber and Haze Correction Computations 

The system shall have the capability to compute the green number 
of each dot for each acquisition and to compute a set of green 
index numbers and a haze parameter for each acquisition. The 
green number for a pixel is the value of the Tas sailed Cap greenness 
for the dot minus the greenness of the soil line. These numbers 
will be used to construct analyst aids in dot labeling. The set 
of green index numbers are ciomulative histogram values of the green 
number with anomalous pixels removed. They will be used to assess 
the condition of the vegetation in a segment on a segment--wide 
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"basis. The haze correction parameter is a measure of the amount 
of haze that was present in the atmosphere over a segment when 
an image of that segment was acquired. These ntraihers will be used 
for R&D purposes. All these computations shotild be performed 
for each acquisition that is loaded into the LACIE Data Base with 
the results output to cards according to the format of 
Section It. 1.4 and the cards delivered to CAMS. Both the green 
ntimber and haze computations use several standard operations 
that will be described first. 


4.1. 4.1 


Recalibrate the Data and Rotate to the Tasselled Cap Coordinate 
System 


Let X = 


where X* ~ relative radiance as provided by 
Landsat MSS band j . 


X4 


4x1 


y 





vhere 




4x1 


relative radiance as provided by 
Landsat-2 MSS band j when processed 
at GSFC with the pre-launch values 
of calibration parameters 



Rj^ and b 


k 


describe the transformation from case k to the Tasselled 


Cap coordinate system where: 


k = 1 denotes the case of Landsat-2 MSS data that was 

processed at GSFC using the pre-launch values of the 
calibration parameters. 

k = 2 denotes the case of Landsat-2 MSS data that was 

processed at GSFC using the post-mission values of the 
calibration parameters . 
k = 3 denotes the case of Landsat-3 MSS data 

^3 ' j J (y - ^ \) 

1=1 
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4.1A,2 


Sun angle correction 

Let 6o = reference solar elevation 

0 = solar elevation of the acquisition "being processed 


Z = 



where Z. 

3 


/sin 0o\ 
\sTh"{J ■/ 


^0 


j = 1, 2, 3, 


The vector Z represents the Landsat data for a pixel that has 
been recalibrated to the Landsat-2 pre-launch calibration, sun 
angle corrected and rotated to the Tasselled Cap coordinate 
system. This is the data that will be used in the screening 
test and in the haze computation. The vector y represents the 
Landsat data for a pixel that has been recalibrated and rotated 
to the Tasselled Cap but not sxm angle corrected. This is the 
data that will be used in the green number computation. 

It . 1 . 4 . 3 Screen out bad data, clouds-, cloud shadows and water, etc . 

The screening process consists of 6 tests; the first 5 are 
used in the haze correction computation and all are used in the 
green number computations . 


Test 1 - A pixel is labeled as garbled data if 


Zi ^ Zi 

4 4, max 


or if Z, < Z, 

4 4 ,min 


or if Z, - 0.093T5 Z, > K 

3 1 3,max 

or if Z^ + 0.18750 Z^ < Z_ 

3 ''^1 3, min 


or if Z„ + 0.1 Z < Z„ 

2 12 ,min 


or if Zp + Z /1.8 > Z 


or if' Z- - Z^/1.2 > Z^ 

2 1 6, max 


Test 2 - A pixel is labeled as cloud if not labeled garbled and 

me 

/ 


Z- ^ Z and Z^ + 0 .1 Z^ ^ Z 
1 cmax 3 1 cmin 


lt-2T 



Test 3 - A pixel is labeled as dxffuse cloud if not labeled 
garbled or cloud and 


2, > 

1 Hjinax 


and 

S * \min 

Test 4 - A pixel is la^beled as vater if not labeled garbled, 
cloud or diffuse cloud and 


and 

and 

and 

and 


Z-, < Z , 

1 wl 


Zg Z^/16 < Z„21 


"i, " "w4 


^ 2 * ^24 


^2 * O'? 2,213!. 


Test 5 - A pixel is labeled as cloud shadow if not labeled 
garbled 5 cloud diffuse cloud or water and 


and 


^2 - \ S - " 's2134 

Zi - 0.4 Z^ < Z^^2 


Test 6 - A pixel is labeled as too bright for use in the soil 
line computation if 


Z-, > Zo 

1 2, max 

4. 1.4. 4 GiTen an acquisition, the green value of the soil line and the 
green index number for the acquisition are computed as follows 

Step 1: Select correct Tasselled Cap transformation and 

Screening Test parameters 

Step 2: Recalibrate data and rotate to the Tasselled Cap 

coordinate system 

Step 3: Apply Sun Angle correction 

Step 4: Screen the data using Tests 1-6 
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step 5: Construct the cumulative histogram for the green 

value without sun angle correction (y^) for all 

pixels that pass the screening test, i.e., 

I 

G(i) Is the nimiher of "good” pixels with 72 > I 

Step 6: Find the green value of the soil line (IS) as the 

smallest yg such that at least m pixels have 
y^ < IS where m = x number of pixels that 

passed the Screening Test. 

NOTE: IS can be negative; however, y^ > -32 for an 

pixels that passed the Screening Test. 

Step T: S is the set {0,1,2,3, ... ,20,22,24, ... ,30,35,^0,45,50} 

Step 8: For JeS, GN(J) = G(lS+J) 

Step 9: GIN(O) - x 100 

Step 10; GIN(k) = x 100 for keS 

Step 11: Bound each GIN(k) to 12 

If GINCk) > 100, set GINCh) = 99 

4. 1.4. 5 Given a pixel with Landsat measurement vector x some 
acguisition for which the green value of the soil line is IS, 
the green number (n) of the pixel is computed as follows; 

Step 1: Select correct Tasselled Cap transformation and 

Screening Test parameters 

Step 2: Eecalibrate Landsat data and rotate to Tasselled 

Cap coordinate system 

Step 3: Compute green nmber for pixel 

n = 72 - IS 

4. 1.4. 6 The haze correction parameter Cy) Is computed as follows: 

Step 1: Select correct Tasselled Cap transformation and 

Screening Test parameters 

Step 2: Recalibrate Landsat data and rotate to Tasselled 

Cap coordinate system 
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step 3: Apply sun angle correction 

Step U: Screen acquisition using Tests 1-5 

Step 5- Compute segment haze diagnostic 


X = 


1 

m 


I 

S 


S = the set of m pixels that passed the Screening 
Test m < 22,932 


Step 6: 


Compute 

k 


“ ' [Ji ""is ] 

” ' [Ji ] ii 


sin 0o 
Sin 0 ' 

sin 6 q 
sin 9 


r V '' AT si 

' iJl ii 


sin 0 


ii. _ V* 


Si- _ Y' 

sin 0 ^ 


If h^ - 2ac > 0, then y = 1000 x ~ - 2ac ) + 300 


a a 


If h - 2ac < 0, then y = 1000 x - + 300 
"" a 


where 


j* = 


n * n 


^1 



* 


“2 

X2 


a.. 


a = 

3 

X3 


_«4 

* 


Lx4_ 



and y* are parameters. 


Jj-.1.5 The system must contain the option to make scatter and 

trajectory plots. The default is to produce both types of plots. 
These plots must he output via line printer as follows: 
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Trajectory plots .- At C&I loading time produce trajectory plots 
of green nxunber vs. time and "brightness vs. time for all 209 dots 
and all acquisitions on the data "base. These plots are to "be 
produced sequentially hy dot number with three dots per plot . 

The brightness plot for a given set of three dots shall be 
directly beneath the green number plot for the same dots and 
there shall be three sets of plots for each page of printout . 

Each page should include segment number and acquisitions on the 
first line. The plots are to show a symbol 1, 2, 3 or * with 
green number (or brightness) read vertically and acquisition date 
horizontally. The symbols 1,2 and 3 refer to the sequential 
order of the three dots represented (l = first dot, 2 = second, 

3 = third) and * records more than one dot with the same green 
number (or brightness) on a given acquisition. The green number 
axis shall range from 0 to 50 with 2 green numbers per line 
(25 lines). The brightness axis shall range from 30 to 100 with 
2 brightnesses per line (35 lines). The horizontal axis shall 
consist of h2 spaces , with each space representing a nine day 
interval starting on September 1 . Consecutive days shall be 
plotted adjacent to one another with the earliest date plotted 
as it falls, and the second day on its right. Each set of 
plots shall be annotated with dot nxmiber of the first dot in 
the sequence, its location (line and pixel) and a sequence of 
characters indicating the labels of the corresponding dots 
(a blank shall be output if no label exists on the data base). 

To allow approximate coordinates to be read off of either plot 
the vertical axes shall be annotated in increments of 10 units , 
and the time axis shall show the approximate start of alternate 
months starting with September. 

it. 1 . 5 . 1.1 Cluster trajectory plots.- Provide the option to produce 

trajectory plots from cluster means of the acquisitions used 
during batch processing. These plots are to be green nvimber vs. 
time and brightness vs. time. They are to be produced sequentially 
by cluster number with three clusters per plot and arranged 
in the same format as the C&I trajectory plots (Ref. it.l.5.l). 

U.l. 5.2 Scatter plots .- Provide the capability to produce scatter plots 
of green number vs . brightness on a fixed scale for each 
acquisition used in batch processing. There shall be two types 
of scatter plots: 

A. One plot per acquisition shall display the category label 
of each labeled dot. Unlabeled dots shall be displayed on 
this plot usxng a separate symbol. Multiple hits shall 
be recorded as the total number of hits rather than as a 
category symbol. 

B. One plot per acquisition shall indicate the category in 

which each dot was classified. In the event that classification 
IS not performed, the category label of the cluster to 
which the pixel is assigned should be used in place of the 
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classifier label. A separate symbol shall be used for dots 
that were thresholded or lie in DO/DU fields. Multiple hits 
shall be recorded as the total number of hits rather than 
as a category symbol. 

The capability to relate dot spectral values to dot grid numbers 
shall be provided as follows. A background xy-grid shall be 
printed on each plot to allow the approximate coordinates of any 
dot to be read off. In addition, the following shall be provided: 
Dot number table - A seq.uential listing by dot grid (not 
random sequence) number giving on the same row: dot grid 
number, dot grid coordinates (scan line n\jmber and sample 
ntonber), analyst label, classifier label and dot green number 
and soil brightness for each acquisition rounded to the 
same precision as used in making the scatter plots. 

Green number table ~ For each acquisition a listing of green 
number, brightness, dot number, grid intersection, analyst 
label, and classifier label. This list shall be ordered by 
increasing green numbers. Dots with the same green number- 
shall be further arranged in order of increasing soil 
brightness. 

The above plots and tables shall be annotated with segment 
' number and acquisition date. 

1^.1. 5. 3 The option shall exist to apply sun angle and/or haze corrections 
to any of the above mentioned plots . The default option shall 
be to produce plots without these corrections. 


vv,‘ 
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SECTION 5.0 

SUBSYSTEM INPUT REQUIREMENTS 


Inputs required By CAMS from the other LACIE' functional elements 
are designated in the following subparagraphs . 

5.1 

5.1.1 

5. 1.1.1 


These images' are to be stored by the ISRRS. Access to approximately 
40 frames/day is required. 

5. 1.1. 2 Lands at MSS segment data .- (Category 1, except as noted; Implementa- 
tion Responsibility: DAPTS) The CAMS performs classification on 

117 line by I96 pixel segments of Landsat data. These segments are 
received from the Data Acquisition Preprocessing and Transmission 
■ Subsystem (DAPTS) via the ISRRS. 

Segments shall be of two types, ordinary segments and designated 
training segments. The acquisition and preprocessing of ordinary 
segments is required in each of fo ur time periods unless such 
acquisition is prevented by cloud cover, technical problems, or 
operational problems. The acquisition and preprocessing .of desig- 
nated training segments is required for each of up to 16 opportu- 
- nities per year unless: 

A. Acquisition is prevented by cloud cover, technical problems, 
or operational problems. 

B. Acquisition occurs within two days after a previous successful 
acquisition. 

The tot^ number of such segments, up to 4800, shall be determined 
by CAS. Of these segments, no more than 960 may be designated 
training segments at any one time; however, the system shall be 
designed such that any segment may be made a designated training 
segment or ordinary segment at any time, as long as this does not 
result in more than 96O active designated training segments. 

Preprocessing of segments shall consist of data quality evaluation, 
calibration, and registration, 
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DAPTS-Landsat Data Requirements 


Regular Landsat imagery .- (Category 1; Implementation Responsibility; 
DAPTS) A 1:1 000 000 scale Landsat (9"3(9") image of each Landsat 
frame used in LACIE is required fo\ir times each growing season. 

This imagery shall be CIR composite of bands 4, 5, and 7. The 
imagery shoiild contain less than 20 percent cloud cover. The 
imagery should be representative of each of the four defined bio- 
logical phases and for each of the 4800 LACIE segments. 
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Data quality evaliiation shall test the data for excessive cloud 
cover, loss of pulse code modulation decommutation synchronization, 
or noise. Data shall be evaluated as good, marginal, or bad. Good 
data will be input to the CAMS. Marginal data will be input to the 
ISRRS for CAMS, and replaced i-f good data are reordered in time. 

Bad data will not be input to the ISRRS. 

Calibration shall adjust all data such that identical spectral 
radiances input to the Landsat MSS result in identical outputs of 
the GSFC preprocessing system. Outputs are regarded as identical 
if the noise introduced by the MSS-GSFC system is negligible in 
comparison with the noise introduced by atmospheric variations and 
other short-term temporal variations in the radiation input to the 
scanner. 

Registration shall align each successive acquisition of a segment 
such that the RFK distance between the ground positions correspond- 
ing to a given image pixel on two separate passes is less than 1 
pixel. Absolute registration accuracy to the surface of the earth 
shall be such that no segment is more than 5 statute miles from its 
nominal position. 

When a sample segment maintains the same nominal position in two 
succeeding crop years, its actual position shall be registered 
Within TBD pixels. (Category 1; Implementation Responsibility: 

DAPTS) . 

Capability to change any JSC defined parameters associated with a 
sample segment, or to order reacquisition of data, shall be provided. 

5. 1.1. 3 PFC gain and bias factors .- (Category 1; Implementation Responsi- 
bility: DAPTS) Gain and bias factors for each segment shall be 

computed in the preprocessing and included as part of the header 
record, for the purpose of controlling the gain and bias settings 
on the production film connector. These factors will be computed 
using an algorithm supplied by CAMS. 

5.1.1.^ Mosiacs .- (Category 1; Implementation Responsibility: DAPTS) 

Black and white uncontrolled 1:1,000,000 scale mosaics of each 
of the LACIE country wheat growing regions are required. These 
mosaics shall be available before the start of analysis in any 
individual country. Each mosaic should be compiled using Band 5 
(the red band) from imagery acquired during the wheat growing 
seasons. Two sets of each mosaic are required, one for the ADP 
country analysis group and one for the image analysis group. These 
shall be on a stable base. 

For each major wheat growing region a north-to-south set of 
uncontrolled strip mosaics is required. This mosaic shall also be 
at a scale of 1:1,000,000, but should be compiled from a CIR com- 
posite of Landsat bands 4, 5, and 7. Each mosaic shall be compiled 
from a single Lands at-? pass, acquired -onc^e during each biological 
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phase using an average or midlatitude of wheat region crop calendar 
to select the actual calendar time to maJse the data acquisition. 

Each image used should contain less than 20 percent cloud cover. 

5.1.2 DAFPS-Eield Data Requirements 

Items A, B, and C of 5»1*2.1 are needed by the end of each calendar 
year for early evaluation of LACIE sites . Item C relates to bio- 
window 1 and data for fall planted crops are needed to coincide 
with FOB processing of Landsat data. These data will also support 
RT&E and accuracy assessment requirements. 

5. 1.2.1 OroTond truth - CAMS .- (Category 1: Implementation Responsibility: 

DAPTS/S&AD-FSO/USDA) 

DAPTS will provide the following data for each LACIE Intensive 
Test Site. ' ‘ 

A. One copy of a 1:24,000 scale rectified B&W photo-mosaic of 
each ITS. The photo-mosaic is to be made from high altitude 
aerial photography acquired once per year, just before wheat 
harvest . 

B. One copy of a 1:24,000 scale transparent overlay registered to 
the photo-mosaic of (a) above and delineating the boundaries 

of each field in the Test Site. Each field will be numbered and 
a tabular listing provided containing the following data; 

1 . Field number 

2 . Field acreage 

3. Crop and variety being grown 

4. Is field irrigated? If so, how? 

5. Is field fertilized? If so, how? 

6 . Planting date 

’These data are to be acquire“d“once per crop year, after spring 
planting and before wheat harvest. 

C. One copy of a 1:24,000 scale field boundary map registered to 
the photo-mosaic of (a) above and delineating the boundaries 
of only those fields containing fall planted crops in each 
Intensive Test Site in Texas, Kansas, Indiana and Washington. 
Each field will be nmbered and a tabular listing provided 
containing the information of item B above. These data are 

to be acquired once per crop year, after fall planting and be- 
fore beginning of dormancy. 

D. One copy of Periodic Observation data acquired at l8 day inter- 
vals from approximately 50 fields in each Intensive Test Site. 
Fields selected for observation must contain both wheat and 
other crop grown in the area. Fields selected must be 40 
acres or larger in size. Adjacent fields planted in the same 
crop and divided by fences or farm roads may be grouped to ob- 
tain a 40 acre or larger area. Data are required from the end 
of dormancy for fall planted crops or spring crop plantings 
‘until harvest. Data are to be acquired concurrent with 
Landsat-2 overpasses • The following data are required in tab- 
ular fo 2 mi. 
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1 . Growth stage 

2. Percent ground cover 

3. Plant height 

4. Surface moisture 

5. Weed growth 

6 . Field operations 

7 . Growth/yield detractants 

8. Stand quality rating 

9. Vertical color photos 

10 . Observation date 

11. General comments on field and site conditions 


5.1.3 


DAFTS-Historical Agricultural Requirements 

A. All of the following data requirements are needed to coincide 
with receipt of any Landsat data at JSC with the following 
exception: Where signature extension is being worked, the data 

are to be provided prior to Landsat acquisitions. The follow- 
ing is a list of countries which require data prior to Landsat 
acquisition: 

Sept 1, 1975 U.S. Great Plains, Winter Wheat 

Oct 1, 1975 USSR, Winter Wheat 

Oct 1, 1975 Canada, Saskatchewan 

Dec 1, 1975 Canada, Remainder 

Jan 1, 1976 U.S. Great Plains, Spring Wheat 

Feb 1, 1976 USSR, Spring Wheat. 


5. 1.3.1 Soil maps .- (Category 1; Implementation Responsibility; DAPTS/USDA) 
Four copies of individual soil maps for the U.S. and each LACIE 
foreign country are required. The areas within each coimtry for 
which soils information is required are those areas used for 
sample segment definition. Maps are to be provided as a positive 
print, at >a scale of 1:5 000 000 or greater (as available), and at 

a known projection. The soil classification system preferred is 
the 7th approximation system, however the genetic (Dokuchev) system 
is also acceptable. In both classification systems there must be a 
sufficient level of detail on the maps provided to allow determina- 
tion of the boundaries of soil suborders . Maps must be reviewed 
yearly and updated as new material becomes available prior to each 
LACIE phase . 

5. 1.3. 2 Topographic maps .- (Category 1; Implementation Responsibility 
DARTS /FSO) 

One copy of topographic maps for all areas containing sample seg- 
ments are required for all LACIE coimtries. In the U.S. these maps 
are required in scales of ±:2k 000 and 1:250 000 for each sample 
segment an,d 1:1 000 000 for areas containing one or more sample 
segments. Substitute maps of 1:62 500 scale are acceptable where 
l:2ll- 000 scale maps are not available. The contour interval shall 
be the standard for the scale map under consideration with the cri- 
terion being to show the greatest amount of detail while maintain- 
ing clarity. 
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For the LACIE foreign countries the maps required will he at, or 
as near as possible to, the scale and detail of the U.S. maps. In 
cases where topographic maps are not available, relief maps showing 
feature changes using hachures, tints or shading will be provided. 

These maps ,(both U.S.. and foreign) will be updated as better in- 
formation becomes available. 

5. 1.3. 3 A priori probabilities .- (Category 3; Implementation Responsibility: 
.DAPTS/USDA) DAPTS has the responsibility of providing raw input data 
for the operation of a computer program to compute the prior proba- 
bility of wheat in a given unit area. The input parameters are TBD. 

A computer program specification, detailing the computational pro- 
cedures to be employed,, the frequency of computations, input data 
sources, output formats and any other required pertinent data and 
information required to develop the probabilities and their associ- 
ated errors will be supplied to DAPTS by CAMS. This program will 
be implemented and operated by DSAD (Required by CAMS; Category 3; 
Implementation Responsibility: DAPTS/DSAD) . 

The various categories of agricultural ground cover to be considered 
are listed in the following subparagraph. 

5. 1.3. 3.1 A priori probabilities : (Category 3; Implementation Responsibility: 

DAPTS/USDA) There is a. requirement to detemine that, if a seg- 
ment of land is cropland, the cover class is a member of one or 
more of the following categories : 

A. Winter wheat - 
‘B. Spring wheat 

C . Wheat 

D. Winter, cereal 

E . Spring cereal 
P. Cereal 

G. Fall plowed crop 

H. Spring plowed crop 

I. Harvested winter wheat 

J. Harvested spring wheat 

K. Harvested spring cereal 

L. Harvested winter cereal 

M. Harvested fall plowed crop 

N. Harvested spring plowed crop 

The estimated error in each of the above probabilities is required. 
These values are to be stored in the electronic ISRRS and are to be 
available to the CAMS whenever the CAMS is operating. An inter- 
active query capability shall- also be provided. 

5. 1.3. 4 Ancillary summaries .- (Category 1; Implementation Responsibility: 
DAPTS/USDA) ' DAPTS shall output to the LPDL for the CAMS one copy 
of an ancillary data sxmmiary for each of the 4800 .sample segments . 
The summaries shall contain information on soils, climate and 
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cropping practices as well as agricultiiral statistical data. The 
suimnaries shall he presented "in the format shown hy the attached 
examples (figs. 5-1 to 5-3). Each summary will be accompanied hy 
its appropriate Locator-Map. Only two general formats are to he 
used - one for the U.S. and one for all foreign countries. The 
Soviet Union summary is an example of the format to he used in 
foreign countries. The text of these formats will he varied to 
meet the potential practices of a specific country. 

The ancillary summaries are required at the county level for all 
counties in the United States containing one or more 5 x 6 n. mi. 
LACIE sample sequent s. For the foreign coxmtries listed below, the 
data are required for political subdivision designated as acreage 
strata for aggregation purposes. 


Country Acreage Strata 


Canada 

U.S.S.R. 

P.R.C. 

India 

Australia 

Argentina 

Brazil 


Census Division 

Oblast 

Province 

District 

Division 

Partido 

Municipality 


The agricultural statistical data required can he determined from 
the ancillary sianmary forms themselves. The cropping practices 
data requires an explanation of the particular information required. 


Cropping practices are defined as any cultural activity that may 
create or alter the external morphological characteristics (appear- 
ance) of a field or crop signature response as detected hy hi^ 
altitude aircraft or spacecraft remote sensors. Examples of crop- 
ping practices are as follows : 

A. Cropping systems or crop rotation practices 

B. Tillage practices 

C. Irrigation practices 

D. Use of fertilizers 

E. Use of special-purpose crops 

1. Catch crops 

2 . Cover crops 

3 . Green-manure crops 
U . Soiling crops 

5. Silage crops 

6. Hiorse crops 

F . Drainage practices 

G. Row directions and widths 

H. Chemical applications (herbicides, defoliants, etc.) 

I . Burning 

J. Field conditions post harvest (stacks, hales, plowed, etc.) 

K. Mowing patterns 

L. Grazing 

M. Harvesting practices 
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-UNITED STATES ANCILLARY SUMMARY 


Coimty (ies). 
Seg, Coord., 


.State. 


.Seg. No. 


.CRD. 


Soil Type/Color. 


Climate (Koppen). 


CROPPING PRACTICES 


Seasonal Precipitation 
Annual Precipitation 


Figure 5-1.- United States ancillary summary. 
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County/State 


Land Use 

Acreage 

(Most Recent Year) 

Percent of Total 
County Area 

A. Total County Area 



B. Total Cropland 



1* Cropland Harvested 



2* Cropland Pastured 



3. All Other Cropland 



C» Woodlands, V»oodland Pasture 



D. All Other Land 

' 



E. Average Field Size (Acres) 

F. Wheat (% Total of County Area) 


II. Conservation Practices (Most Recent Year) 

A, Irrigated Land (% County Area) 

B. Contour (Acres) Strip (Acres) Terrace (Acres) 


III. Crops and Agricultural Lands (Absolute Acreage and i of Total County Area) 


Crops 

(Include All) 

Most Recent Year 

Second Most Recent Year 1 

Acreage 



Acreage 

Percent 

Percent 

Irrigation 

^Vheat 




■ 




Figure 5-1 Concluded, 
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SEGMENT NO 


SOVIET UNION 
ANCILLARY SUMMARY 

S . S . R . OBLAST/KRAY 

SEGMENT COORDINATES LAT. LONG. 

MAP NO . - - • SCALE 

SOIL TYPE/COLOR CLIMATE (KOPPEN) 

WMO NO. 

I. CROPPING PRACTICES ' 

Figure 5-2.- Soviet Union ancillary stunmary. 
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Land Use 

A. SSR 

B. Oblast/Kray 

1. Cropped Area 

a. Net area sown 

b. Area sown more than once 

c. Irrigated area 

2. Fallowed Area 

.a. Fallowed less than I'year 

b. Fallowed 1-5 years 

c. Fallowed more than 5 years 

3. Non-cultivated Area (a+b+c) 

a. Forest 

b. Permanent pasture/grazing 

c. Other (marsh/desert/ 

urban, etc.) 


Area in 1000 Ha. Percent as 

Most Recent Year Indicated 




















STATE SEGMENT No. 

CRD & No. DATE 


CROP CALENDARS 

PERCENT OF AREA IN DEVELOPMENT STAGE BY SPECIFIED DATE 


CROP %ar 
100 

WINTER 75 

WHEAT 50 

25 

75 

EA 


^ARC 

:h 

APRIL 


MAY 




JUN 


JULY 

AUGUST 

SEPTEMBER 

OCTOBER 
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75 

50 

25 

75 

SO 

25 

75 

SO 

25 

75 

50 

25 

75 

50 

25 

75 

50 

25 

75 

50 

25 

75 

50 

25 

75 

50 

25 

75 

50 

25 

75 

50 

25 


• YEARS OF DATA - AVERAGE 

(NO ADJUSTMENTS FOR WEATHER MADE) 

• SOURCE OF DATA 

• ESTIMATED ACCURACY OF DATES AS THEY WOULD APPLY TO ANY GIVEN YEAR 

± DAYS 


Figure 5-3.- Crop calendars. 
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The following is a further amplification of the types of informa- 
tion required. 

A. Farm and field patterns: Data concerning the farm and field 

patterns in each wheat region will he required. 'In addition, 
representative schematics will he needed to illustrate "typical 
field patterns and land uses in the wheat regions. The follow- 
ing schematic shows the land use on a farm. This illustrates 
the' type of schematic required. 



B. Farm sizes: Information on the size of the farm and the indi- 

vidual fields (hy crop type) is required, hut a narrative ex- 
plaining the data is hi^ly desirable to the Image Interpreter. 
An example is this narrative excerpted from the Yearbook of 
Agriculture, 1958. 

"The average commercial wheat farm in West-Central Kansas has about 
TI 5 .acres of land. Of this, 525 acres are cropland and 19 O acres 
are native pasture. About 330 acres of the cropland are seeded to 
wheat. As about two-thirds of it is seeded on summer- fallowed land 
the average farm has about l40 acres of cultivated summer-fallow 
each year. 

An average of about 40 acres, or one-sixth of the seeded wheat, 
fails to make a crop because of winter killing, drought or other 
hazard. If it rains after the wheat has failed some of this 
abandoned wheat land is planted to grain sorghum the following 
spring. The possibility of planting a catch crop of sorghum on 
abandoned wheat land is a unique characteristic of wheat-growing 
in the Southern Plains . . " 

On the average, about 70 acres of grain sorghimi and 4o acres of 
forage sorghimi are planted each year. . . 

The average farmer has about 40 head of cattle and no other live- 
stock on a commercial scale." 

C. Crop rotation systems: Generalized information on the crop 

rotation systems will be needed. Three examples of this type 
of data are as follows : 

1. Crop rotation systems in the Corn Belt: 

a. 3-year rotation of com, oats, red clover timothy sod 

b. 4-year rotation of corn, soybeans, oats and clover 

c. 5-year rotation of two years of corn, one year of oats, 
and two years of alfalfa-bromegrass . 


2. One-year cropping systems: 

a. Winter wheat followed "by annual lespedeza 

b. Winter rye grass ^alternated with a summer row crop such 

‘ as corn . 

3. Cropping practice in Yangtze Valley, China: 

a. Cotton is hand broadcast in the wheat field one month 

t 

prior to the wheat harvest. 

Storage of these data on a sample segment basis in the LPDL is re- 
quired. Approximately 20 such sets of data per day will be 
accessed. 

5. 1.3. 5 , Historical Crop Calendars .- (Category I; Implementation Eesponsi- 

bility: DIPTS/USDA) DAPTS is required to provide a mean historical 

crop calendar for each sample segment to assist the Analyst Inter- 
preters (A.i.'s) in image interpretation. Each calendar shall de- 
scribe the progress ,of wheat and all other crops which, taken to- 
gether, constitute at least 95 percent of the cultivated area for 
which information is available. The year or years of historical 
data used in the development of the calendar shall be stated. Data 
sources shall also be stated. For foreign countries the calendars 
shall be prepared for areas similar in size to a U.S. Crop Report- 
ing District (CRD). . In the .U.S, data will be required at the CRD 
level . 

The calendars for each crop will be presented as percent area in a 
specified growth -stage at a -specified date. An example of the de- 
sired data product output format is attached. 

5. 1.3. 6 Acquisition windows .- (Category 1; Implementation Responsibility: 
DAPTS) The CAMS is responsible for determining the actual Landsat 
MSS data acquisition phases based upon estimated best acquisition 
phases plus Landsat- MSS data already acquired. 

The designated acquisition windows are defined by the CAS in terms 
of crop growth stage. These windows then input to the Yield Esti- 
mation Subsystem (YES) . Estimates of dates delineating the growth 
stages on the basis of crop growth models , historical meteorological 
data, and/or other information will be used. These estimates are 
transmitted, -via the ISRRS, to the DAPTS for action and to the CAMS 
for conciurence . ,If the CAMS fails to concur, it must provide 
acquisition windows to replace those provided by the YES . 

To accomplish the above , the CAMS requires the dates delineating 
all crop acquisition windows for all segments as transmitted from 
YES to DAPTS. The ISRRS will maintain accessible records of all 
such transmissions . 
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5.1.4 YES Requirements 

5. 1.4.1 Crop calendar updates .- (Category 1; Implementation Responsibility: 
YES) A crop calendar is required for each acreage strata for which 
a sample segment exists. The following defines a crop calendar as 
required to support the image interpretation. 

A ’’point plot curve" (from DAPTS) providing the expected time peri- 
ods of the seedbed preparation, planting, phonological stages, and 
harvest of the significant crop types within a specified region for 
each biological year that LACIE is in operation. YES will provide 
a tabular update to these data every two weeks during the growing 
season . 

For each field condition or biological stage the earliest date, the 
latest date, and that period in time when approximately 50 percent 
of a given crop has reached a defined stage should be included. 

(Ref. 5. 1.3. 5 for format.) 

This type of calendar is required for wheat, its confusion crops 
(small grains), and other crops such that calendars are provided for 
95 percent of the ciiltivated area in the neighborhood of the sample 
segment . 

Updates to these calendars shall be provided to indicate early or 
late attainment of various growth stages of wheat and other crops 
in response to real-time meteorological data. The format shall be 
the same as used for historic calendars when more crops than wheat 
are, updated. 

Storage of these data on a sample segment basis in the ISRRS is 
required. Approximately 20 such sets of data per day will be 
accessed. 

Individual crop calendars should be provided for each acreage 
strata for which a segment is designated for CAMS AI analysis . 

Sample segment numbers with coiintry and locality descriptors should 
be included in the title data for each calendar. 

Source data for the calendar construction shoTold be clearly stated 
on each calendar, including the data of the latest meteorological 
adjustments applied. 

Interpolated data should be distinguished from USDA or other 
statistical sources. 

5. 1.4. 2 Meteorological summary .- (Required by CAMS; Category 1; Implementa- 
tion Responsibility: YES) Meteorological data will also be utilized 

to provide weekly updates of the weather being experienced in speci- 
fic regions . Such data will be similar to that being provided for 
the United States in the Weekly Weather and Crqp Bulletin; i.e., 
normal moisture and temperature , wetter and warmer than normal , 
drier and warmer than nomal, and drier and cooler than nonnal. 
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5.1.5 


These be provided in tabulations which will apply to the 

sample segments. A narrative, which would flag areas where signa- 
ture extension might not be successful due to some parameter other 
than temperature and/or precipitation, will also accompany the 
above. Such other parameters include wind damage, hail damage, or 
other anaomalies occurring within the yield strata. 


CAS Requirements 

The following information is reg.uired a^: least one month prior to 
the start of any LACIE operational phase : 


5 .1.5.1 Sample allocation . - ■ (Category 1 ; Implementation Responsibility: CAS) 
The nominal locations of all sample sequent s shall be indicated on 
transparent overlays suitable for .registration upon USAF Operational 
navigational Charts (OEfC's) or Landsat uncontrolled mosaics. Sample 
segment center coordinates shall be latitude and longitude in degrees 
and minutes. These overlays shall also indicate the achieved Landsat 
ground tracks , LACIE frame boundaries , and LACIE overlap zones . 

Three copies of each such overlay are required. 


5-1. 5. 2 


Landsat acquisition- rates (Category 1; Implementation Responsi- 
bility: CAS) Histograms indicating' the expected acquisition rates 

in segments per month are required separately for winter and spring 
wheat for each co\mtry. 


ALSO for each LACIE designated country provide individual segment 
numbers for each predicted acquisition date during the wheat pheno- 
logical development stage. 


5. 1.5. 3 


Sample segment histogram .- (Category 1; Implementation Responsibil- 
ity: CAS) Histograms indicating the distribution of number of 

sample segments per unit area for the spring and winter wheat 
regions of each coimtry are required. 


5.1.5.^ Flagged area output .- (Category 1; Implementation Responsibility: 
CAS), Sample segment classification that results in .more than 
20 'percent change in the previous biophase percentage wheat area 
will be returned to CAMS for classification verification. 


5.1.6 


ISRRS Requirements 


5. 1.6.1 Initial imagery .- (Category 1; Implementation Responsibility: 

ISRRS ) Initial imagery products shall consist of color composite 
and black and white imagery generated from Landsat segment data. 
This imagery shall be produced, at the nearest practicable scale, 
to 1:125 000, such that constant spacing between pixels occurs in 
the crosstrack direction and constant spacing between scan lines 
occurs in the along-track direction. 


Black and white images shall display the value of each pixel in one 
channel on one acquisition as a set of levels of gray ranging from 
black to white. j The .number of shades of gray, n, shall be 8, 16 , 
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or 256, as a user option. Data values shall be assigned to gray- 
levels as follows : 

1st gray level (black) corresponds to 

y = 0to(^)-l 
2nd gray level corresponds to 

3rd gray level corresponds to 

5^ = ^ 1 3 . 1 

last gray level (white) corresponds to 

where 


y = A (X + B) + C 

where X is the pixel data value for the channel of interest, A 
is a scale factor, B is bias 1, and C is bias 2. A, B, and C 
may be used selected; default -values are indicated below. 

The user shall have the option to select any channel on any acqui- 
sition date for imagery generation; default values are indicated 
below. 

Color images shall display the value of each pixel in a selected 
band in each of three colors (red, green, and blue). Each color 
shall be available in a range of brightness ranging from black to 
saturation of the film. The mmiber of color levels, n, shall be 
the same for each color; the user shall have the option of select- 
ing between n = 8, 16, or 256. Data values shall be assigned to 
each color level as follows : 

1st color level (black) corresponds to 

y = 0 to y = - 1 

2nd color level corresponds to 



last color level (satiaration of film) corresponds to 


vhere 


y = A (X + B) + C 

and X is the pixel data value in the selected channel, A is a 
scaling factor, B is the second "bias factor, and C is the over- 
all bias. 

A, B, and C may be user selected. Default values are indicated 
below. 

The user shall have the option to assign any single channel/ 
acquisition date combination or pair thereof to any color in the 
image . Default values are indicated below. 

All Landsat MSS imagery shall be generated at the same scale. 

All default imagery products shall be produced within 2k hours of 
JSC receipt of the associated segment data. User req.uested imagery 
shall be produced within 2k hours of user request. 

Flexibility to change default values \mder configuration control 
shall exist. 

Gain and bias values for all Landsat MSS segments are supplied in 
the header record for all tapes received from GSFC on a per channel 
bias. These values shall be used as indicated below and shall be 
stored for future use . 

All products shall have a superimposed grid at every 10th line and 
pixel starting at line one and pixel one. The grid shall not ob- 
scure more than one-half the area of any pixel. 

The above capabilities are operational. 

5. 1.6. 2 Landsat MSS tape products .- (Category 1; Implementation Responsi- 
bility: ISRRS) The products listed b.elow are operational: 

5. 1.6. 2.1 Screening film (Product Number 7)t Screening film shall be pro- 
duced for every acquisition of every segment. It shall be a mod- 
erate gain simulated color IR product using MSS bands U, 5> and 7 
and a monochrome moderate gain product for band 6. 

MSS band h shall be assigned the color blue in the color composite; 
MSS band 5 shall be assigned the color green; MSS band 7 shall be 
assigned the color red; and MSS band 6 shall be used to create a 
monochrome image . 
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In each ease, the overall bias, C, iised for a band shall be taken 
from the tape header and multiplied by a TBD constant. Bias B 
shall be zero. The gain. A, shall be the value taken from the 
tape header multiplied -by the above TBD constant. 

5. 1.6. 2. 2 Color composite positive (Product Number l): This product shall be 

produced for the first acguisition in each biological phase of all 
currently designated training segments. It is a high gain simu- 
lated color IR product using Landsat MSS bands it , 5 j and J . 

The color blue shall be assigned to f©S band it. The bias B shall 
be 0. The bias C shall be taken from the GSFC header record. The 
gain shall be taken from the GSFC header record. 

The color green shall be assigned to MSS band 5* Gains and biases 
are as above. 

The color red shall be assigned to JKS band J. Gains and biases 
are as above. 

5 .1.6. 2. 3 Color composite positive /negative (Product Number 2): This product 

shall be produced for the same acquisitions as the color composite 
positive'. 

The color red shall be assigned to MSS band 5* ^he bias B shall 
be zero. Bias C shall be taken from the GSFC header. Gain 
shall be taken from the GSFC header. 

The color blue shall be assigned to MSS band 6. The bias B shall 
be zero.’ Bias C shall be the bias value in the GSFC header. Gain 
shall be the gain value in the GSFC header. 

The color green shall be assigned to MSS band J. The gain and bias 
shall be as immediately above. 

5. 1.6. 2. U Black and white images (Product Nmber 4); A black and white image 
shall be produced for each band of the first acquisition in each 
biological phase of each currently designated training segment. 

These will be hi^ gain images . 

The biases B shall be zero. The biases C shall be obtained from the 
GSFC header. The gains shall be obtained from the GSFC header. 

5. 1.6. 2. 5 Temporal color composite: This product shall be produced upon 

request in each biological phase, if an earlier phase has been 
acquired. It is a hi^ gain product allowing the interpreter to 
see temporal effects in a single image. 

Eleven different acquisition sequences contain two or more acquisi- 
tions. ' Color and gain assignments for these sequences are TBD. 
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5. 1.6. 2. 6 Image header; Each LACIE Landsat MSS image shovild he accompanied 
hy the following header information on the film. 

Job ID 

System ID 

Tape library 

Date tape generated 

Sensor 

Orbit 

Product requested 

Sample segment ID, File number 

Channels requested 

Bias factors B 

Landsat scene 

PFC run date 

PFC tape ID 

Output mode 

Humber of pixel hits 

Scaling factors A 

Date of sensor data acquisition 

PFC color balance factors 

Color codes 

Polarity 

Repeat of pixel 

Repeat of scan 

5. 1.6. 3 ADP imagery .- (Category 1; Implementation Responsibility: ISRRS) 
(See section 6.1.4) 

5. 1.6. 4 Optional Landsat imagery .- (Category 3; Implementation Responsi- 
bility: ISRRS) TBD 

5.1.Y Research Requirements 

(Cat. 1; Impl. Resp. : Research) TBD 

5.1.8 Test and Evaluation Requirements 

(Cat. 1; Impl. Resp.: Test and Evaluation) TBD 
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SECTION 6.0 

C.AMS OUTPUT REQUIREMENTS 


Outputs from CAMS ■which •will he provided to the other LACIE func- 
tional elements or "which are output hy the CAMS and reinput to the 
CAMS are defined in the following subparagraphs. 

FULL-UP SYSTEMS 

Inputs to the ADP System 

(Req.*d by CAMS; Cat. 1; Impl. Resp. CAMS/NASA-DSAD ) 

Semiautomatic mode control inputs .- The Classification Applications 
Program semiautomatic mode requires control input. Inputs to the 
processing control data base shou3.d be by card. The analyst will 
process data in the semautomatic mode by accomplishing the follow- 
ing functions -via card input. 

A. Analyst name and run ID. 

B. Recognition segment numbers. 

C. Training segment numbers. 

D. Channel maps between recognition segments and training segments. 

E. Channels to be used for clustering (default is all). 

F. Choice of one of two sets of clustering default parameters 
and option to select chaining for each cluster nm. 

G. Each cluster nm initialized by input classes and subclasses 
with analyst defining the final class name for each clustering 
run. 

H. Classes or sxibclasses to be used for each Bhattacharyya chain- 
ing run. 

I. Normal computation of statistics by class option. 

J. Utilization of external statistics from checkpoint tape. 

K. Generation and -writing of statistics on output tape. 

L. Application of sun angle correction and mean level adjustment 
for various desired recognition-training segment pairs for sub- 
class statistics. Including the option to classify the training 
segment with unaltered statistics . 

M. Classes to be ■used in channel selection. Option to generate 
film. Option to distinguish between interclass distances and 
intraclass distances. Number of channels to be selected. 

N. Option to set category thresholds and a priori probabilities 
for this job. Both classification and feature selection. 

O. Definition of categories to be classified to class rather than 
category level . 

P. Run purpose. 

Q. A listing of the categories to be bias corrected and the name 
(if any) of the category to be excluded from bias correction. 

R. Choice of whether or not to produce spectral plots and choice 

of type of trajectory plot to be produced. Defa-ults are to 
make spectral plots using all acquisitions in the Data Base; in 
the case of consecutive day acquisitions take the first alcqui- 
sition. The default type of trajectory plot is option 1 (ref. 
SEC. It. 1.5). ' 



S. Choice of whether or not to employ feature selection. The 

default shall be to by-pass feat\ire selection unless a special 
input is made by the user. 


6 . 1 . 1 . 2 , 


Updates to fields data base .- The CAMS requires the capability to 
update the fields data base for each segment which contains the 
field definitions, a priori probabilities, and category thresh- 
old values for that segnent. This is accomplished by means of the 
following interim products; 

A. Control cards containing field definitions, ^ £rioi^ probabili- 
ties, and/or threshold values (interim output product). 

B. Tape containing field definitions, a priori values, and/or 
threshold values . This tape is generated by the interactive 
mode of the Classification Applications Program (interim out- 
put product ) . 

C. Multifile tape containing field definitions. This tape is 
created by an offline program and unloaded from the data base. 
The segment information in the header of the tape is 
specified by the user or defaulted to the current segment 
identification . 


6. 1.1. 3 


Updates to dot data base .- The dot data base can be updated by 
means of the following products . 

A. Control cards containing dot number, type, and label. 

B. Multifile tape containing dot mmiber, type, and label. This 
tape is created by an offline program and unloaded from the 
data base. The segnent information in the header of the tape 
is specified by the user or defaulted to the current 
segment identification. 


6.1.2 Data Status Reports 

(Req's by CAMS; Cat. 1; Impl. Resp.: CAMS/RASA-DSAD) 

6. 1.2.1 Indexed data report .- An indexed data report shall exist for every 
sample segment. This report shall be updated following any opera- 
tion resulting in a change in, deletion of, or addition of param- 
eters in the report or upon query. 


This report shall indicate the following parameters for each 
sample segment : 

1. Sample segment niimber 

2 . Country 

3. Spring or winter wheat sample 
li-. Region, zone stratum 

5. Latitude and longitude of nominal center point 

6. GSFC registration and data quality parameters (reference scene, 
orbit, Landsat scene, percent cloud cover, registration- 
sharpness, P/B ratio, sun angle, sync-loss -lines, storage 
status, processing date) 

7. List of acquisition dates for segment 

8. Designated training segment or ordinary segment 

9. Ciirrent definition of biological windows 

10. History of updates to GSFC data base 
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6 . 1 . 2. 2 

6.1.3 
6. 1.3.1 


11. Accession information for all data received from GSFC for this 
segment 

12. Latest run ID. 

Acquisition status matrices .- To accomplish manual training seg- 
ment selection, access to the following items will he required: 

1. Acquisition status at JSC of all segments'. ' 

2. Acquisition status at GSFC of all segments. 

3. nominal acquisition dates for each segment. 

Classification Besults 


Classification summary report .- (Required hy CAS, CAMS; Category 1; 
Implementation Responsibility: CAMS/NASAyDSAD ) The processing of 

each segment hy the classification processor in either the semi- 
automatic or interactive mode shall include a classification sum- 
mary report . This report shall include the following : 

A. 1. Date and number of RFC classification map generated for 

this processing segment. 

B. 1. Threshold values used and actual number of pixels 

thresholded/ class , category and subclass . 

C. 1. A priori probabilities used per category. 

D. Classes and subclasses used on data control cards and for 
processing' functions. Sources of all subclasses, i.e., cluster- 
ing, B-chaining, etc. 

S. (Transformed) band means and variances for each class or sub- 
class . 

F. Classified percentages for each category C^, Cg,..,, C to be 
bias corrected where : ^ 

C = name of category to be bias corrected, 1 < q < p. 


p = number of categories to be bias corrected. 

= n umb er of pixels _ in DU fields. 

= number of pixels thresholded in classification 
= number of pixels classified as Category A 
B = 22,932 - 

WpQ = number of pixels in all DO fields 

A = the name. of a category whose pixels are to be excluded 

from classification percentages. A shall be a user input 
with default A=X. 

I \ 

X “ usually an unident if ia'ble category (clouds/cloud ^shadows). 


= the number of pixels classified in category ' 

N - - the ninnher of pixels classified into a category other 
^ than the categories 





-tl - 

= :^x 100; %C = ““ x 100; %I)0 = x 100. 

q B ’ p +1 B B 

Wvmiber of -pixels in DO fields 

196*117 - mmiber of pixels in category 
tmidentifiaBle - nuinber of pixels in 
DU fields 


2. Corrected percentages for each category C^, C2,...jC^. 
For an index, r, 1 2 r ^ p+1 



where 


m = mmher of bias correction dots labeled C by the 
analyst and by the classifier 


m = ntimber of bias correction dots labeled C by the 

<1>P+1 q 

analyst .and as some category other than A, C^, 

C by the classifier. 

2’ ’ p 

m = Wvimber of bias correction dots labeled C by the 
^ classifier. ^ 


m , = number of bias correction dots that were labeled by 
^ the classifier as some category other than A, C^, 

^2 5 • • • J P • 

The summary shall be output such that a sequential listing 
of type 1 dots is output first, followed by all dots (other 
than type l) classified C^. 

Note, any bias correction dots that were thresholded in 
classification or that fall in DU areas shall be excluded 
from m 

q^r 


a 




m /m 
q»r' r 


0 


m _/m 
q,p+l' p+1 


if one or more bias correction dots 
were classified C for q=l,2..., p and 
r=l , 2 , . . . , p 

if no bias correction dots were classi- 
fied C i . e . m =0 
r r 

for q=l,2,..., p and r=p+l 
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A dot processing summary shall he provided for each segment 
processed (whether classification is performed or not) . 

In general, this summary sjiall contain the following infor- 
mation sequentially list'ed'-'with respect to dot nxmibers: 

Dot number 

Dot pixel, row number v 

Dot classification category (as determined hy classifier if 
classified) 

Dot cluster label (label of cluster to which dot was 
assigned by clustering algorithm) 

Dot category label (if specified by analyst) 

Dot type (S-starting vector,. L-labeling vector, B-bias cor- 
rection, DO, DU) 

Dot green number 
Dot bri^tness 

The summary shall be output such that a sequential listing 
of starting and labeling dots is output first, followed by 
all dots (other than starting or labeling' dbts ) classified 
C- then Cp sequentially through C , followed by dots 

classified other than, _A, C.. , Cp, — , C followed by DU, 
threshold, then DO. ^ > 

In addition, the following, se^ent sunm^y information shall 
be output to the dot summary report. 

P(C ) through. P(C ) 
i P 

through 

throu^ 

“ 2.1 “ 2 . P+1 

throu^ 


a T throu^ a . ^ , 
P,1 P>P+1 


For each q 5 -p, a p+1 dimensional vector V whose 
is ^ 



component 


V 




”n 

2 

a ( 


‘ — X 100 

m - 1 
r 
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where Z = I. 6 U 5 . 


Number of pixels in DU fields 

196*117 


I. So^^rce of training data: 

1. Training segpient(s), date(s) of acquisition. 

2 . Training fields . 

3 . External statistics. 

J. Channels used. 

K. Channel map used to relate training segment channels to recog- 
nition segment channels . 

L-. Identity of recognition segment (s). 

i!. Acquisition dates of recognition and training segments, includ- 
ing all flaggings of marginal data. 

N. Whether run in semiautomatic or interactive. 

O. Playback of analyst semiautomatic mode inputs. 

P. Playback of analyst interactive mode inputs. 

Q. Identifying run manber, date of run. This item shall appear on 
all individual products. 

R. Biological phase. 

S. A classification summary similar to the ERIPS classification 
siimmary, generalized to include subclasses, for all test and 
training fields in recognition and training segment (s), and 
for the entire recognition training segment, giving pixels 
classified into each subclass. 

T. As item S, except by class. 

U. As item S, except by category. 


The following reports shall be included with the classification re- 
port when applicable: 

V. A fields definition report similar to that generated by the 
ERIPS, but including subclass name for each field. 

¥. A clustering report including channels used for clustering, 
cluster popvilation, class matches, cluster names, a detailed 
clustering report, interdistance clustering report, intradis- 
tance clustering report, cluster parameters, definition of 
cluster maps generated for PFC and also output on microfiche. 
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X. 


A chaining report including fields /clusters chained, an intra- 
distance report (B-distance; initial and final), initial number 
of subclasses and subclass identities, final number of sub- 
classes and subclass identities, number of pixels per chain. 

Y. A featxire selection report including channels available, chan- 
nels selected, interclass separability, total separability with 
all channels, total separability with selected subset of chan- 
nels, summary of film output of PFC. 

Z , User name . 

AA. For interactive runs, if interactive field changes are made, a 
complete definition of' fields, a priori probabilities, and 
threshold values used for classification shall be reported. A 
field boundary overlay shall also be generated and referenced. 
BB. Diagnostics and incorrect inputs. 

CC . The user shall have the capabilities to print selected pattern 
recognition reports. 

6. 1.3. 2 Electronic outputs .- Electronics outputs are as follows: 

6.1. 3.2.1 Ready access electronic outputs for sample segments: Items 1 to 

30, defined in section U. 1.1. 8. 4, classification output products, 
as being saved on magnetic tape for each classification run shall 
be stored in the ISRRS physical library for ready access use by 
CAS, 


6. 1.3. 2. 2 Spectral value unload: The capability shall be provided to unload 

spectral values for all 209 grid intersection dots for any subset 
of* the acquisitions on the data base to either tape or cards. In- 
puts shall be segment number, acquisition date(s), and output mode. 
Multiple segments may be output to one tape. If no acquisition 
date is input all acquisition dates are tinloaded. Default output 
mode is tape. In tape output mode, the spectral information for 
the dots shall be output as an image in Universal format. In card 
■ output mode , the spectral information for the dots shall be output 

in a standard format along with the following parameters from the 
header record: segment number, acquisition date, solar elevation 

and the green value for the soil line in each acquisition. 

6.1.4 ADP Imagery 

(Req'd by CAMS; Cat. 1; Impl. Resp. : CAMS/ISRRS-DSAD) 

Imagery products shall consist of transparencies, at the approximate 
scale of 1:125 000. The exact scale shall be nominally identical 
to that of other LACIE applications imagery products. ADP image 
outputs shall be overlaid with a line and pixel ground identical to 
that used for initial imagery. 

6. 1.4.1 Field verification overlay .- (Product 'Rtmiber 12).-' (Required by 

CAMS; Category 1; Implementation Responsibility: CAMS/ISRRS-DSAD) 

This shall be an image indicating the bovuidaries and identities of 
all fields in a given segment . It shall be produced for a segment 
following any occasion on which the fields data base for that seg- 
ment is updated, or following any interactive inin in which field 


6-7 



^ definitions are changed prior to classification. This requirement 
is operational. 

6.1,k.2 Cluster map (Product Number 9) •- (Required by CAMS; Category 1; 
Implementation Re spons ibility : CAMS / ISRRS-DSAD ) . 

Color coded conditional and unconditional cluster maps (ref. sec. 

4. 1.1.1}-) shall be output in image form each time a segment is batch 
processed. The boundaries and identities of all field(s) shall be 
indicated by a technique which does not obscure any pixels. It is 
permissible to utilize multiple sheets to accomplish this. This 
cluster map shall be produced for all clustering runs. This map 
should also contain an identifying number. This requirement is 
operational. 

6. 1.4,3 Classification map (Product Ntimber 8).- (Required by CAMS; 

Category 1; Implementation Responsibility: CAMS /ISRRS-DSAD ) The 

classification map shall indicate wheat vs harvested wheat vs other 
vs unidentifiable vs thresholded. The boundaries and identities of 
all field(s) shall be indicated by a technique which does not ob- 
scure any pixels . It is permissible to utilize multiple sheets to 
accomplish this. The classification shall be produced for all 
classification jobs. This map should also contain an identifying 
run number. This requirement is operational. 

The grey shades or colors used in the classification maps shall be 
consistent by category between segments, yet retaining contrast. 

For example, wheat on Segment A should be represented the same as 
wheat on Segment B. CAMS shall have the capability to specify 
default values. 


6. 1.4. 4 Dot Overlay of Analyst Labels .- The overlay shall indicate the cur- 
rent label of each analyst labeled dot and shall be at the same scale 
as the LACIE AI Imagery. It shall be produced for a segment follow- 
ing any occasion on which the dot data base is updated. 

The segment summary information shall be output to the Dot Summary 
report whenever cluster bias correction is performed (paragraph 
4.1. 1. 4. L) . 

6.1.5 Outputs to DAPTS 

(Req'd by DAPTS; Cat. 1; Impl. Resp. : CAMS/LACIE ASVB) 

6. 1.5.1 Updates to GSFC data base orders (Computer Card Coding Sheet s).- 
CAMS shall implement the requests for changing the data acquisi- 
tion and preprocessing order to GSFC. These changes will include 
specifying a need for additional processing of a sample segment for 
a given biological phase, or deleting sample segments, changing 
film processing flags, and changing color code information. The 
change will be specified by making the appropriate entries on a 
computer card coding sheet, using the format contained in the JSC 
card format documentation. These changes will then be submitted 
through the EOD Data Manager. 
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6 , 1 . 
6 . 1 . 

6.1 


6.1 

6.1 

6.1 


5.2 List of training segments .- CAMS shall specify, by segment nvraiber, 
those segments designated as training segments. This designation 
shall occur for up to 9^0 of the sample segments . 

5.3 FFC gain and Mas computation algorithm .- CAMS shall provide the 
algorithm to he implemented by GSFC in computing the gain and bias 
settings for the PFC. 

6 Data Base Transaction Query Report 

(Req’d by CAMS; Cat. 1; Impl. Resp. : CA^K/NASA-DSAD) A data b^se 

transaction report shall be generated for all alterations of stored 
field definitions, a priori probabilities, threshold values. This 
report may also be produced upon request to determine data base 
status. If the job produces a change in data base status, a copy 
will also be output for archival purposes. The report shall indi- 
cate the following: 

A. R\m number 

B . Date 

C . Transactor name 

D. Data base initial status (field names, classes, types coordi- 
nates ; ^ priori probabilities ; threshold values ; date of most 
recent update). 

E. Data base final status (as above) indicating all changes in 
status . 

F. Accession information for field boundary overlay if applicable. 

G. Accession information for all previous transaction reports. 

The total nvunber of reports shall not exceed 100,000. 

T Resource Parameter Reports 

T.l Manual resources parameter reports .- (Required by SPE; Category 2: 
Implementation Responsibility: CAMS/NASA-ASVB) Throughput and 

resources information will be collected for SPE or an as required 
basis. 

.8 Spectral Plots 

Each time a segment is batch processed spectral plots (ref. sec. 
4.1.5) shall be output to a line printer if the spectral plot option 
has been chosen. 
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SECTION T.O 

INTERFACE REQmREMENTS 


Not applicable 
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SECTION 8.0 


OPERATIONAL REQUIREMENTS AFFECTING THE SYSTEM DESIGN 

The following requirements apply to all organizations implementing 
systems or procedures in support of CAMS. 

8.1 THROUGHPUT REQUIREMENTS (Required Ly ’CAS.;. Cat., l) . ^ 

The system -shall he capable of delivei^ing to CAS the following 
nvimber of satisfactorily processed segments after allowance is made 
for rework, accuracy assessment, quality control, data base main- 
tenance,, research, test, evaluation, and software development: 

1 April 1975 '7 -31 'December 1975 100 sample, < 

— segraents/week - . - 

1 January 1976 - 30 September 1976 170 sample 

— segments/week 

1 October 1976 - 31 March 1976 200 sample 

— segments /week 

8.2 RESPONSE REQUIREMENT 
(Req’d by CAMS; Cat. l) 

The turnaround time for all products associated with a classifica- 
tion processing operation shall be not more than 2h hours. 

8.3 RELIABILITY REQUIREMENTS 
(Req’d by CAMS, Cat. l) 

' The system reliability and throughput shall be such that the back 

of CAMS inputs or outputs accumulated due to failure of any system 
component is less than one day^s peak processing capability* 

8.1i SECURITY REQUIREMENTS 

A LACIE Security Plan shall be prepared by each organization 
designated with implementation responsibility. The plan for each 
implementing organization shall define the specific measures that 
will be utilized by that organization to comply with the LACIE 
seciority requirements. 

The LACIE security requirements shall be defined by the USDA and 
will be forwarded to the implementing organizations upon receipt 
by the LACIE Project Manager. 

Each implementing organization shall submit their plan for approval 
to the LACIE Level III Change Board within 90 days after their re- 
ceipt of the USDA requirements. 

8.5 DELIVERY OF PRODUCTS PAGE BLANK NOT 

(Req'd by CAMS; Cat. l) 

Delivery of all CAMS inputs and outputs to place of use is required. 
For those products to be analyzed by EOD personnel, the place of 
use is building 17 » JSC. 
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8.6 QUALITY ASSURANCE REQUIREMENTS 

(Reel’d "by CAMS; Cat. l) 

A LACIE Quality Assurance Plan is teing prepared under the direc- 
tion of the Quality Asstirance Manager. 

8.T OTHER OPERATIONAL REQUIREMENTS 

(Req’d hy CAMS; Cat. l) 

The system shall be so designed that any segment can be processed 
on any day, regardless of the processing status of other segments. 

The system shall be designed to minimize the resources expended in 
preparing ADP inputs. Specifically, CAMS analysts shall not be 
responsible for card piinching operations. 


8-2 



SECTION 9.0 

CAMS VERIFICATION REQUIREMENTS 


9.1 

9.1.1 

9.1.2 


9.2 

9.3 

9.3.1 

9.3.2 


DATA SYSTEMS VEEIFICATION 

(Req'd by CAMS; Cat. 1; Impl. Resp. : CA&B/NASA-DSAD) 

Design Reviews 

Preliminary and critical design reviews of the CAMS data system 
shall be conducted. The CAMS manager shall have the authority to 
submit Review Item Dispositions against portions of the design 
affecting the CAMS. 

Verification Testing 

The verification of the CAMS data system shall be accomplished 
jointly by the EOD and the DSAD according to the following outline: 

A. The implementing organization shall prepare a test plan accept- 
able to the DSAD and to the CAMS manager. 

B. Following acceptance of the test plan, the implementing 
organization shall prepare a test procedure acceptable to the 
DSAD and to the CAMS manager. 

C. Following acceptance of the test procedure, the implementing 
organization shall perform the test, witness, ed by the subsystem 
manager or his designate(s) . A document of acceptable test 
completion shall then be prepared. Whenever practicable, accep- 
tance tests will utilize data sets previously analyzed in the 
test and evaluation program. 

D. Quality assurance will approve on test plans and test proced- 
ures 

ANALYST CERTIFICATION 

(Req'd by CAMS; Cat. 1: Impl. Response: CAMS/LACIE-ASVB) 

All analysts performing classification processing are required to 
be certified in accordance with the LACIE CAMS Analyst Certifica- 
tion Procedures Document, LACIE-00300. 

PROCEDURES VERIFICATION 

(Req'd by CAMS: Cat. 2; Impl. Resp.: CAMS/NASA-ASVB) 

Operational Procedures 

All operational CAMS procedtires shall be documented and verified 
against their intended objectives, as stated in the Procedures 
Purpose. This verification shall be accomplished according to a 
plan and procedure acceptable to all organizations affected by the 
procedure to be verified. 

Analysis Procedures 

All CAMS analysis procedures shall be documented and tested before 
implementation. The documentation shall consist of a functional 
procedure which depicts the analysis in detail flow charts and a 
detailed procedure which technically describes the functional flow 
charts. 
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SECTION 10.0 
BESEABCH REQUIBEMENTS 


(Required "by CAMS; Cat. 4; Impl. Resp. Research) 

TBD after development of LACIE Research, Test, and Evaluation 
Plan. 


blank Not 
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SECTION 11.0 
TEST MP EVALUATION 

(Required "by CAMS; Cat. U; Impl. Resp. : Test and Evaluation) 

TED after development of LACIE Research, Test, and Evaluation 
Plan. 
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APPENDIX A 


REQUIREMENTS FOR 
CAMS INTERFACE REPORT 

AS PRODUCED FROM THE CAMS /CAS INTERFACE TAPE 


A sample of the required report and an explanation of the data 
in each part of the report are attached. 

Three options are required for the production of this report. 
They are as follows : 

A. The option to produce the entire report for .all 
segments on a given tape, 

B. The option to produce the entire report for specific 
segments on a given tape, and • 

C. . The option to produce reports for all segments on a 

given tape, hut for each segment, produce only the 
reports in Articles I, IJT^III and IV of the 
attached exp>lanation and sample unless one of the 
computed PCC's for that segment is less than 80, and 
produce the entire report for that segment otherwise. 


BtAWK fi&T 


A-1 




EXPLANATION OF THE REPORT 


I. HEADER INFORMATION . This includes 

A. DPAR nuiaber 

B. Record ID 

C . Segment number 

D. A list of the acq.uisitions used in processing 

E. The tape number from which the report was generated. 


II. CLASSIFICATION AND CLASSIFICATION BIAS CORRECTION RESULTS . 
This includes 

A. Type 1 Dot Label/Classification Matrix 

This matrix shows the user assigned label and classified 
category of each starting and labeling dot. Restilts for 
each dot are displayed according to dot line and pixel 
nvimber. The user label is to be printed first, followed 
by a slash and then the classified label. Dots that are 
thresholded in -classification or that lie in DO or DU 
fields should have their classified categories reported 
as TH, DO, or DU accordingly. 

B. Type 2 Dot Label/Classification Matrix 

This matrix shows the user label and classified category 
of each bias correction dot and the classified catetory 
only for unlabeled dots. This information is to be 
reported in the same format as the Type 1 Label/ 
Classification Matrxx (Article II B). 

C. Bias Correction Classification Report 
This report lists: 

1. For each category of interest 

a. Category classified popiilation 

b. Category classified percentage 

c . Category bias corrected percentage 


A-2 



d. Category variance 


e. Category random sample estimate computed by 


RSE- ( ) =‘^he number of bias correction 

dots labeled in the specified 
category and classified in any 
category of interest (not 
thresholded) 




the total number of bias 
correction dots which were 
i classified in a category of 
interest without thresholding; 


xlOOx / 22932-population DO-population unidentifiable 
V 22932-population unidentifiable 


2. For DO, threshold and unidentifiable categories 

a. total population of the category in question 

b. the percentage of the category in question on 

the segment. . . 


3. If ¥ and S are both categories of interest, then 
compute and report items in Article II C.l above 
for grain category ,as' follows i 

Population of grain category is the siim of the 
populations of category W and category S 

’ • ( ' 
b. Classified Percentage is the sim of the classified 

percentages 'Of category ¥ and category S 
c- Corrected Percentage is 



+ (l-“lT) 



X 100 


ot 

where G = The number of bias correction dots 

labeled either ¥ or S and classified 
(not -thresholded) either ¥ or S 
The total number of bias correction 
dots labeled with any categoiy of 
interest and classified (not 
thresholded either ¥ or S, 
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= The number of bias correction dots 
labeled in any category of interest 
other than W or S and classified 
(not thresholded) in any category 

of interest other than W or S 

The total number of bias correction 
dots labeled with any category of 
interest and classified (not thres- 
holded) in any category of interest 
other than W or S, 

N is the population of category W, 

is the population of category Sj 

is the sum of the populations of all 
^ categories of interest other than ¥ or S, and 

Base is 22932 less the population of unidentifiable. 

Variance is 


N +N 
w s 


Base 


X IOC 


acCl-Og) 


+ 



Base 


X 100 


I13-1 


¥here n is the nimiber of bias correction dots 
^ classified (not thresholded) in 
category ¥, 

n is the number of bias correction dots 
^ classified (not thresholded) in 
category S, and 

n is the number of bias correction dots 
^ classified (not thresholded) in a 

category of interest other than ¥ or S. 



e. Random Sample Estimate is 


n^+iis+nj ^ 


N +N +N, 
w s 5 

Base 


xlOO 


Where is the number of bias correction dots 
labeled category ¥ and classified in 
any category of interest and (not 
thresholded) , and 

g is the number of bias correction dots 
labeled category S and. classified 
in any category of interest (not 
thresholded) . 

If any denominator is zero in these com- 
putations, the corresponding term should 
be zero. 


D. Alpha Value Matrix . 

The rows and columns of this matrix represent categories 
of interest. Each entry reports the Alpha value for the 
category of the row of the entry vs . the category of the 
coliann of the entry. 


E. 


PCC ' s . These are computed-as - follows : 


PCCl = 


i 


\ 


Number of starting and labeling dots both 
labeled and classified (thresholded dots not' 
included) in the same category of interest. 
Number of starting and labeling dots which 
were labeled and classified (thresholded 
dots not included) in any of the categories 
of interest . 


X 100 


2 . 


PCC2 


The ninnber of bias correction dots both 
labeled and classified (excluding thres- 
holded dots) in the same category of 

interest . 

The number of bias correction dots 
(excluding those thresholded) which were 
labeled and classified in any of the 
^ categories of interest. / 


\ 




X 100 


I 


\ 


If both categories ¥ and *S are categories of interest 
then compute PCC ' s for the grain category as : 
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3. PCCGl = 


The number of starting and labeling dots 
labeled W or S and classified (not thres- 
holded) either ¥ or S + the number of 
starting and labeling dots labeled in a 
category of interest other than W or S 
and classified (not thresholded^) in a 
category of interest other than ¥ or S 
The number of starting and labeling dots 
which were labeled and classified (not 
thresholded) in a category of interest. 

PCCG2 is computed the same way as PCCGl except that 
bias correction dots are used rather than 
starting and labeling dots. 

F. Bias Correction Vectors (Labeled; Classified) . 

This report lists for each ordered pair, (A, B), of 
categories of interest, the number of bias correction 
dots labeled in the first category (here category A) 
and classified in the second category (here category B) . 
This number is reported as A;B. The report also 
includes the sum of all these entries as total. 


G. The values Dw and Ds are computed as follows: 


8 8 


8 


Dw = + ^p=i ^ category ¥ not used'- 

8 8 8 

Ds = ^p=i ^q=i UpiSpSq. + ^p=p VpSp or 0 if category S not used 


8 8 8 8 8 8 

where {Rpq}^ {Upq.} ^ , and {Vp} _ are 

p-1, q.-l, p-1, p-1, q-1, p-1 

constants (TBD) . 


= Bias corrected estimate for category ¥ 
= Bias corrected estimate for category S 
Wg = Machine estimate for category ¥ 

Sg = Machine estimate for category S 
= Random estimate for category ¥ 

= Random estimate for category S 


100 
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= Variance of bias corrected estimate of W 
= Variance of bias corrected estimate of S 

W = W;W 

5 W:S + S:S + F:S 

_ S^S 

W:S + S:S + N:S 

Wg = Sg = PCCl 

y = S = 

T T W:K + S:W + U:N 
Wg = Sg = PCC2 

The ntunbers W:¥, ¥:S, etc. are as described in article II F 
above. If category W or category S is not a category of 
interest then each entry involving that category is assumed 
to be zero. 

III. CLUSTERING AM) CLUSTER BIAS CORRECTION REPORT 


The format for this report is the same as the format for the 
bias correction report {article~Il) . All entries are made using 
the category of cluster to which dots are assigned rather than 
classified label. Dots' -which are unas signed to a cluster are 
treated as thresholded dots. 

IV. SEPARABILITY REPORT 

This is a listing of feature selection Battacharyya separabil- 
ities for all available acquisition combinations . 

V. CLASSIFICATION DATA 


This reports the cluster name (subclass name), threshold value, 
and normalized apriori value for each cluster (subclass) used 
to train the classifier. It also lists the channels used in 
classification. 
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VI. 


CLUSTER REPORT 


VII. 


VIII. 


This report shows the following: 

A. Clusters Generated : The total ntmiber of clusters generated. 

B. Options Used : A listing of options used in clustering. 

Specifically sun angle normalization, DO/DU exclusion, or 
min/max radiance exclusion. 

C. For each of the clusters generated, a listing of cluster 
name, labeling dot, the clustering distance from the 
labeling dot to the cluster mean, the user assigned labels 
of dots assigned to that cluster, the cluster brightness 
and green number. This is followed by listing of the dots 
assigned to the cluster, their user assigned label and 
type and their clustering distance to the cluster mean. 

This list should be ordered by ascending clustering 
distance. 

STATISTICS REPORT 


This report lists, for every subclass used to train the 
classifier, the subclass name, population, and mean and 
standard deviation in each channel. 

DOT REPORT 


This reports for each dot in the order present on the tape of 
the following information: 

A. Dot Number 

B. Line 

C . Pixel 

D . Dot type 

E. User Label 

F. Label of cluster to which dot was assigned 

G. Classification category (Th if thresholded) 

H. Brightness / For each acquisition \ 

I. Green Number \ used in classification/ 
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